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RESUMO 


MOTA, Júlia de Araújo. Estudo da teoria de lubrificação com parametriza- 
ção diferenciada da geometria e aplicações em mancais hidrodinâmicos. 
2020. [137]. Dissertação (Mestrado em Informática) - PPGI, Instituto de Matemá- 
tica, Instituto Tércio Pacitti de Aplicações e Pesquisas Computacionais, Universi- 
dade Federal do Rio de Janeiro, Rio de Janeiro, 2020. 





As equações de Reynolds são amplamente usadas para entender a geração de 
pressão na teoria dos lubrificantes. Entretanto, avanços tecnológicos causaram di- 
ferenças nos parâmetros operacionais e no material utilizado, justificando a neces- 
sidade de uma revisão nas simplificações tradicionais. Por esse motivo, o presente 
trabalho tem o propósito de revisitar as simplificações feitas nas equações de Navier- 
Stokes para o escoamento entre dois cilindros com uma pequena espessura de filme 
de óleo lubrificante. Através de uma análise adimensional, os termos dessas equações 
serão mapeados e ordenados por importância no fenômeno físico para esse tipo de 
situação. Uma parametrização eficaz da geometria será proposta, viabilizando uma 
descrição mais detalhada do problema e sua adaptação para outros contextos. No 
limite das simplificações, encontra-se uma equação diferencial parcial elíptica, resol- 
vida pelo método de diferenças finitas centradas e que tem como solução o campo 
de pressão entre os cilindros. Para ilustrar o sucesso da abordagem, o modelo será 
aplicado em mancais hidrodinâmicos (componentes de máquina que fazem interface 
entre as partes móveis e fixas de um sistema rotativo), sendo calculados o campo de 
pressão e alguns parâmetros decorrentes do mesmo, como por exemplo, coeficientes 
de rigidez e amortecimento, para os quais será proposto um método diferente do 
tradicionalmente utilizado pela literatura. Com as respostas obtidas, serão realiza- 
das análises sobre o impacto da variação dos principais parâmetros geométricos e 
operacionais do problema. A partir das facilidades oferecidas pela parametrização 
da geometria, apresentam-se outras duas configurações de mancais hidrodinâmicos: 
mancais elípticos e com desgaste, avaliando suas respectivas respostas e promo- 
vendo um estudo comparativo entre as geometrias. A modelagem exposta neste 
texto, bem como todas as suas simulações foram desenvolvidas para integrar Ross- 
Rotordynamics, uma biblioteca em Python para análise rotodinâmica, disponível na 
plataforma GitHub. 


Palavras-chave: Teoria de Lubrificação, Parametrização da Geometria, Simulação 
Numérica, Mancais Hidrodinâmicos, Coeficientes Dinâmicos. 


ABSTRACT 


MOTA, Júlia de Araújo. Lubrication theory study with differentiated geo- 
metry parameterization and applications in hydrodynamic bearings. 2020. 
137]p. Dissertation (Master in Computer Science) - PPGI, Institute of Mathema- 
tics, Tércio Pacitti Institute of Computational Applications and Research, Federal 
University of Rio de Janeiro, Rio de Janeiro, 2020. 


Reynolds” equations are widely used to understand pressure generation in lubri- 
cant theory. However, technological advances have caused differences in operational 
parameters and in the material used, justifying the need for a review of traditional 
simplifications. For this reason, the present work aims to revisit the simplifications 
made in the Navier-Stokes equations for the flow between two cylinders with a small 
thickness of lubricating oil film. Through a dimensionless analysis, the terms of these 
equations will be mapped and ordered by importance in the physical phenomenon 
for this type of situation. An effective parameterization of the geometry will be 
proposed, enabling a more detailed description of the problem and its adaptation to 
other contexts. At the limit of simplifications, there is a partial elliptical differential 
equation, solved by the centered finite difference method and whose solution is the 
pressure field between the cylinders. To illustrate the success of the approach, the 
model will be applied to hydrodynamic bearings (machine components that inter- 
face between the moving and fixed parts of a rotating system), being calculated the 
pressure field and some parameters resulting from it, such as stiflness and damp- 
ing coefficients, for which a method different from the one traditionally used in the 
literature is presented. With the answers obtained, analyzes will be carried out on 
the impact of the variation of the main geometric and operational parameters of the 
problem. Based on the facilities offered by the parameterization of geometry, two 
other configurations of hydrodynamic bearings are presented: elliptical and worn 
bearings, evaluating their respective responses and promoting a comparative study 
between the geometries. The modeling exposed in this text, as well as all its sim- 
ulations were developed to integrate Ross-Rotordynamics, a library in Python for 
rotodynamic analysis, available on the GitHub platform. 














Keywords: Lubrication Theory; Geometry parameterization; Numerical Simula- 
tion; Hydrodynamic Bearings; Dynamic Coefficients. 
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1 INTRODUÇÃO 


Estudar o escoamento de um fluido é fundamental para o avanço em diversas 
áreas da engenharia, fato que ganha relevância quando relacionado a componentes 


de máquinas rotativas. Isso ocorre devido à grande utilização de tais equipamen- 





tos no mundo contemporâneo, presentes em contextos antes não imaginados. Um 
grande exemplo no cenário atual é sua utilização na exploração do pré-sal para ex- 


tração do petróleo, submetendo estas máquinas a condições adversas por conta dos 





grandes gradientes de pressão e variações de temperatura, além de diversos tipos de 


camadas a serem perfuradas. Quanto mais adversas são as condições físicas, mais 





importante é analisar o escoamento do fluido, possibilitando entender sua influência 


no funcionamento da máquina. 


Apresentada no final do século XIX, a equação diferencial amplamente utili- 
zada na atualidade para a teoria dos lubrificantes foi desenvolvida por 
(1886). Essa equação descreve a distribuição de pressão do filme fluido no espaço 
anular de um rolamento e é aplicada a diferentes conjuntos específicos de circuns- 
tâncias, gerando diversos tipos de equação de Reynolds. Segundo DOWSON| (1962), 
essa multiplicidade pode mascarar a influência de fenômenos físicos associados à ge- 
ração de pressão. O autor, já naquela época, aponta diferenças entre teoria e prática 
no que diz respeito a parâmetros operacionais e materiais utilizados. Atualmente as 
condições operacionais são ainda mais discrepantes. Isto significa que simplificações 
tradicionais deixaram de ser condizentes com a prática, sugerindo uma nova inves- 
tigação minuciosa nos modelos usados para descrever a velocidade e a pressão nos 


componentes do rotor. 


Aperfeiçoar os modelos supracitados significa ter a chance de obter mais 
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informações sobre o campo de pressão e, consequentemente, de sua influência nas 
vibrações, campo de velocidade e rotação para que o sistema mantenha estabilidade 
dentro de padrões aceitáveis. Por esta razão, alguns componentes merecem ter seus 


modelos matemáticos revisitados, como no caso de mancais. 


Fluido Lubrificante Eixo 


“ 






Mancal 


Figura 1.1: Esquema de um mancal hidrodinâmico. 


Fonte: (RAMOS, Pois) 


CASTRO (2007) define mancais como elementos de máquina que fazem uma 


interface entre as partes móveis e fixas em um sistema rotativo. Se um mancal é 
hidrodinâmico, possui um fluido lubrificante entre as peças rígidas cuja finalidade 
é, de acordo com [ALVES (2011), substituir o atrito seco entre as partes pelo atrito 
viscoso, reduzindo a temperatura de funcionamento, o atrito e o desgaste das super- 
fícies. Por viabilizar altas cargas em altas velocidades, os mancais hidrodinâmicos 
são largamente utilizados em turbomáquinas. Além disso, possui longa durabilidade 


operacional uma vez que não existe contato entre as peças. 


LF 


1.1 Revisão da Literatura 


A maioria dos trabalhos relacionados à teoria da lubrificação são datados do 





final do século XIX, em função da Revolução Industrial e o consequente desenvol- 





vimento das máquinas. Os mancais também tiveram origem neste mesmo período. 
Alguns pesquisadores da época empenhavam-se em resolver o problema do atrito en- 
tre os eixos e suportes das máquinas, o que era causa de grandes perdas energéticas e 


elevados níveis de calor. Este problema os estimulou a buscar resoluções através da 


lubrificação, desenvolvendo métodos teóricos ou experimentais. (CASTRO, 2007). 


De acordo com SAN ANDRES! (2010a), a lubrificação por um filme fluido é 


um fenômeno hidrodinâmico caracterizado por um lubrificante fluindo no espaço es- 
treito entre duas superfícies. Os pioneiros nesse estudo foram os britânicos Osborne 
Reynolds (1842- 1912) e Beauchamp Tower (1845-1904), tal qual o russo Nicolai P. 
Petrov (1836-1920), cujos trabalhos continuam válidos nos estudos atuais. (MEN- 


DES por 


De acordo com |PINKUS! (1987) em sua revisão histórica da teoria da lu- 


brificação hidrodinâmica, estes três homens, em poucos anos e independentes uns 





dos outros, descobriram e formularam o mecanismo da lubrificação hidrodinâmica 


e estabeleceram suas bases como um ramo da ciência da engenharia. Estes pesqui- 





sadores perceberem em comum o processo da lubrificação como resultado de um 
filme fluido que separa duas superfícies sólidas, e não em decorrência da interação 


mecânica entre elas. Este é o aspecto fundamental da lubrificação hidrodinâmica. 


A elaboração de conceitos iniciais sobre a lubrificação teve início em 1883, 
quando Petrov formulou a relação entre a força de atrito e os parâmetros de um 


mancal. Tower, em sua segunda publicação Second Report on Friction Experiments 





(1885), expôs as distribuições de pressão de medidas também em um mancal. Mesmo 
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sendo os responsáveis pelo conceito de lubrificação hidrodinâmica, Petrov e Tower 
obtiveram seus resultados apenas de maneira experimental. O desenvolvimento da 
base teórica foi elaborado por Reynolds, em 1886, com a publicação On the Theory 
of Lubrication and Its Application to Mr. Beauchamp Tower, na qual explicou os re- 


sultados obtidos por Tower e acabou estabelecendo a teoria moderna da lubrificação 
hidrodinâmica. (MENDES, 2014). 


y 





Figura 1.2: Mancal deslizante. 


Conforme indicado por/ORON: DAVIS; BANKOFFI| (1997), a teoria da lubri- 


ficação pode ser ilustrada de maneira mais simples quando considerado um mancal 
de deslizamento lubrificado por um fluido viscoso, como indicado na Figura (1.2). 
Nessas peças de máquina o fluido viscoso é forçado a entrar em um canal convergente, 
criando forças de pressão verticais que podem ser usadas para suportar grandes car- 
gas e, portanto, reduzir o desgaste. Segundo o autor, matematicamente, a teoria da 
lubrificação pode ser vista como uma exploração da disparidade entre duas escalas 
de comprimento: a espessura característica do filme h e uma escala característica do 
comprimento L. O principal requisito para a teoria da lubrificação é que a proporção 
e= h/L seja muito pequena, isto é, e < 1. Desta maneira, Reynolds realizou sim- 
plificações da equação de Navier-Stokes, estabelecendo a equação diferencial para o 
perfil das pressões que atuam entre as duas superfícies. Tal equação ficou conhecida 


como Equação Reynolds e é amplamente empregada até os dias atuais. Segundo 
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PINKUS| (1987), apresenta o seguinte aspecto: 


o [30p) O (500) oh 
Es U E) e 9 [h 3 == Gu (Us + Ui) 9x +av |, E) 


onde x é a coordenada na direção de U, z é a coordenada normal a U (direção do 








comprimento), A é a espessura de filme de óleo, p é a pressão hidrodinâmica, u é a 
viscosidade do fluido, Up e U, as velocidades lineares do mancal e do eixo rotacional, 


respectivamente, e V é a velocidade radial. 


De acordo com CASTRO! (2007), o trabalho de Reynolds ainda introduziu 


muitos outros conceitos novos para os pesquisadores da época, como por exemplo, 
a definição de folga radial e a relação da pressão com o fenômeno de cavitação 


nas partes divergentes dos mancais. Por ser uma equação diferencial parcial não 





homogênea com coeficientes variáveis, a equação de Reynolds tornou-se um grande 
desafio para os pesquisadores da época. Tratava-se de uma equação difícil de resolver 
analiticamente e, mesmo quando resolvida para casos específicos, os resultados eram 
complicados de usar. O autor ainda afirma que uma das limitações para a resolução 
dessa equação foi a falta de informações acerca das condições de contorno necessárias 


para sua Integração. 


A partir desse momento, na tentativa de obter alguma solução, vários auto- 
res passaram a aplicar esta equação em um conjunto específico de circunstâncias, 
desconsiderando alguns fatores físicos ou até mesmo restringindo a casos específicos. 
O próprio Reynolds estabeleceu a aproximação para mancais infinitamente longos, 
para os quais considerou a variação de pressão na direção x muito maior que na di- 
reção z, desprezando a segunda parcela de sua equação (1.1). Contudo, não obteve 


sucesso uma vez que a solução era aceita apenas para mancais de carregamento leve. 


(PINKUS) [1987). 


Nesta área, destaca-se 5LOMMERFELD! (1904), que obteve uma expressão 
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analítica explícita com as simplificações para mancais infinitamente longos. No 


entanto, as condições de contorno utilizadas pelo autor desconsideravam algumas 


condições físicas, destoando da realidade. Segundo MENDES)! (2014), condições de 
contorno mais realistas foram apresentadas por [SWIFT (1932), e em |STIEBER 


(1933), sendo utilizadas até os dias atuais para o cálculo de mancais com carga 


constante. 


Em oposição à ideia da aproximação para mancais infinitamente longos, a 


aproximação para mancais infinitamente curtos ocorreu pela primeira vez em IME 
CHELL| (1929), quando o autor sugeriu negligenciar o primeiro em vez do segundo 
termo na equação de Reynolds (1.1). Contudo, apenas em[OCVIRK! (1952), foi apre- 


sentada uma solução detalhada e completa para o problema. Nela a relação L/D 
(comprimento por diâmetro do mancal) é válida para valores menores ou iguais a 


0.9, que é a faixa de design da maioria dos mancais modernos (PINKUS, 11987). 


Esta solução é amplamente utilizada até os dias atuais. 


PINKUS! (1987) indica que antes da chegada dos computadores, [AMERON; 
WOOD (1949) eSASSENFELD; WALTER! (1954) apresentaram soluções numéricas 


para mancais finitos. Ambos os métodos envolveram uma quantidade extremamente 
grande de cálculos, sendo considerados, por suas maneiras de resolver a equação 


diferencial, os precursores da nova era caracterizada pela utilização de computadores. 


Os computadores eletrônicos de alta velocidade ocasionaram uma verdadeira 
revolução no que diz respeito à tentativa de encontrar uma solução adequada para a 
equação de Reynolds (1.1). foi o primeiro a utilizar computadores 
modernos neste intuito e com as condições de contorno adequadas, apresentando 
soluções para mancais circulares, elípticos e de três lóbulos. A característica de maior 
destaque deste trabalho foi a realização de soluções para problemas de diferentes 


geometrias, além do fato adicional de principiar o uso de soluções numéricas para 


2 
determinar as forças hidrodinâmicas em mancais. (PINKUS) |1987). 


Como resposta aos primeiros trabalhos utilizando soluções computacionais, 
em pouco tempo surgiram muitos modelos e métodos de resolução a fim de solucionar 
o problema de lubrificação hidrodinâmica. Um método aproximado para a solução 
da equação de Reynolds foi desenvolvido por (1976). 
(1981) aplicou Diferenças Finitas à equação, propondo a utilização de vários métodos 
iterativos para a resolução do problema, como os de Jacobi e Gauss-Seidel. Já em 


CAVALCA; CATTARUZZI (2001), os autores sugeriram soluções para diferentes 


tipos de mancais também por Diferenças Finitas. 


Várias análises do escoamento de fluidos com diferentes simplificações são 
propostas na literatura, por vezes para aperfeiçoar a teoria sobre o fenômeno da lu- 


brificação, mas também na busca por maiores informações em casos experimentais. 


É interessante ressaltar os estudos de |DIETZEN; NORDMANNI (1987) e, trinta anos 
depois, de/WAGNER. et al.| (2017) que, na busca pelos coeficientes rotodinâmicos, 


notaram a necessidade de antes obter aproximações para o escoamento e a distribui- 
ção de pressão. O primeiro trabalho resolve por Diferenças Finitas um escoamento 
turbulento a partir equações de Navier-Stokes associadas a um modelo de turbu- 
lência, ao passo que o segundo se baseia numa simplificação, utilizando medidas 


experimentais. 


Segundo |PINKUS! (1987), os coeficientes rotodinâmicos se tornaram um ele- 


mento básico nos estudos de mancais desde que Stodola, em 1925, constatou que um 
mancal não é um suporte rígido, mas sim um conjunto de molas e travas, cujas ca- 


racterísticas influenciam os aspectos críticos do rotor e o comportamento dinâmico. 





Tais coeficientes são amplamente utilizados em sistemas rotativos sustentados por 
mancais hidrodinâmicos. Podem ser de rigidez (kra, kry, Kyr € kyy), que relacionam 


a força do filme fluido ao deslocamento do eixo, ou de amortecimento (Cry, Cry, Cyr € 


Pa 


Cyy), relacionando as forças com a velocidade de vibração. [LUND (1987) promoveu 


uma revisão acerca do conceito desses coeficientes dinâmicos, propondo métodos 
para calculá-los. Entretanto, grande parte da literatura sobre o tema desenvolve os 
coeficientes dinâmicos considerando alguma das aproximações citadas anteriormente 


para a equação de Reynolds, como |FRISWELL et al. (2010), que apresenta equações 


para os coeficientes de rigidez e amortecimento em um mancal curto. 


1.1.1 Estudos relevantes 


Algumas obras foram fundamentais para a construção deste trabalho, ser- 


vindo de inspiração e fundamentação teórica. No que diz respeito à teoria da lubri- 


ficação são relevantes os estudos feitos por |PINA; CARVALHO! (2006), a dissertação 


delANDRADE/) (2008) e o projeto de graduação deiIQUEIROZ (2018). PINA; CAR- 
VALHO! (2006) buscam analisar o escoamento associado a perfuração de poços de 


petróleo e gás, considerando uma excentricidade e validando seu modelo através 
da comparação com outro já conhecido. (2008), por sua vez, utiliza a 
teoria da lubrificação para modelar o escoamento monofásico no interior de bom- 
bas de cavidades progressivas (BCP), reduzindo o custo computacional sem perder 
a legitimidade dos resultados. Já analisou o escoamento de um 
fluido entre dois cilindros, obtendo resultados para o campo de pressão, velocidade 


e vazão. 


FOX; MCDONALD; PRITCHARD/! (2014) de forma simples e didática, esta- 


belecem toda a fundamentação teórica acerca da teoria da lubrificação, enfatizando 
os conceitos físicos da mecânica dos fluidos e os métodos de análise que se iniciam a 
partir dos princípios básicos. Por essa razão alicerçou o modelo que será apresentado 


neste texto. 
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Outro trabalho de extrema importância trata-se da conhecida coletânea de 





notas de aula registradas pelo professor San Andrés, especialmente as notas de 1 à 


5. Em sua primeira nota, LAN ANDRES! (2010a) desenvolve, a partir das equações 


fundamentais da mecânica dos fluidos, as equações clássicas para o movimento de 


um lubrificante em uma película fina. Na nota 2, apre- 
senta a derivação da equação de Reynolds na teoria clássica da lubrificação. A 
nota 3, (2010c), introduz a análise do escoamento de fluido em um 
mancal cilíndrico, enquanto a 4, (2010d), apresenta os mancais de 
comprimento longo e curto. Por fim, em sua quinta nota, 


realiza uma análise no equilíbrio estático do mancal, define os coeficientes dinâmicos 





e aponta fórmulas para as aproximações no mancal curto. 


A dissertação de mestrado de |IRAMOS)! (2019) permite uma melhor compre- 


ensão do fenêmeno da cavitação. O autor avalia os efeitos da conservação de massa 
no modelo de cavitação do mancal hidrodinâmico através de uma modelagem da 


equação de Reynolds pelo Método dos Elementos Finitos. 


Os trabalhos de |FRÊNE et al. (1990), HAMROCKI (1991) e |ISHIDA; YA- 
MAMOTO! (2012) destacam-se, para este texto, no que é relativo às aproximações 


da equação de Reynolds. Os três autores apresentam conceitos básicos da geometria 


de mancais hidrodinâmicos e utilizam a equação aproximada para mancais curtos 


e infinitamente longos. |FRÊNE et al.| (1990) ainda realizam análises acerca da ra- 
zão L/D (comprimento por diâmetro do mancal), eISHIDA; YAMAMOTO| (2012) 


discorrem sobre as forças do filme fluido e rotores elásticos. 


Por fim, a dissertação de mestrado de IMACHADO! (2011) e a pesquisa de 
MACHADO; CAVALCA! (2015) apresentam soluções para mancais de diferentes 
geometrias. IMACHADO! (2011) analisa características dinâmicas e operacionais 


de diferentes configurações de mancais radiais hidrodinâmicos: um cilíndrico, um 
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elíptico e um trilobular. Já IMACHADO; CAVALCA! (2015) expõem um modelo 


numérico que caracteriza o mancal com desgaste em sua parede, analisando sua 


influência na resposta dinâmica do rotor. 


1.2 Objetivos 





Tendo em vista os avanços tecnológicos e as consequentes modificações sofri- 
das pelas máquinas rotativas, estas passaram a operar com parâmetros diferentes de 
quando a teoria da lubrificação foi desenvolvida. Em outras palavras, as tradicionais 
simplificações da literatura deixaram de ser condizentes com a prática, havendo a 
necessidade de revisitar conceitos no intuito de contemplar novos efeitos físicos an- 
tes negligenciados. Além disso, grande parte dos trabalhos atuais retrata avanços 
baseados nas conjecturas da teoria da lubrificação adaptadas para casos específicos. 
Com isso, a tarefa de encontrar registros detalhados de como essas simplificações são 
feitas torna-se não trivial, dificultando ainda mais o processo de inclusão de novos 


parâmetros. 


Por esta razão, este trabalho tem como objetivo geral revisitar a equação 


de Navier-Stokes, avaliando os efeitos físicos para o escoamento de uma pequena 





espessura de filme de óleo lubrificante entre dois cilindros. A finalidade é promover 
um mapeamento dos termos das equações, ordenando-os por relevância no fenômeno 
físico para esse tipo de situação. Esse processo possibilitará investigações futuras, a 
fim de reincorporar termos, antes negligenciados, e avaliar seus respectivos impactos 


na compreensão do fenômeno. 


Outro objetivo é viabilizar uma descrição do problema mais detalhada e dife- 





rente de como normalmente é feita na literatura sobre o tema, de forma a desenvolver 


uma parametrização da geometria que permita adaptar facilmente o modelo a ou- 
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tros contextos. Intenciona-se possibilitar a inclusão de variações na geometria ou a 
adaptação para diferentes componentes de máquinas rotativas apenas pela descrição 


do problema, sem que haja a necessidade de alterar a modelagem proposta. 





Para promover um primeiro estudo, o presente texto tem como uma de suas 
metas específicas encontrar numericamente o campo de pressão, de modo a reduzir o 
custo computacional envolvido no processo de resolução das equações que descrevem 


este campo. 


Com a finalidade de ilustrar o sucesso da abordagem, o campo de pressão será 
calculado para mancais hidrodinâmicos. O propósito é comparar os resultados deste 
modelo com o mais utilizado pela literatura (aproximação para o mancal curto), 
obtendo numericamente forças, posição de equilíbrio e coeficientes dinâmicos, cal- 


culados tipicamente através do negligenciamento de termos do modelo simplificado. 


Vale ressaltar que as simulações que serão utilizadas neste trabalho foram 
feitas para integrar Ross-Rotordynamic4l] uma biblioteca de código aberto escrita 
em Python para análise rotodinâmica, por 2020), no qual a autora 
deste texto possui coautoria. Em Ross, o presente modelo corresponde à parte de- 
nominada Fluid Flow, responsável pela simulação do filme de fluido em mancais 
hidrodinâmicos que retorna ao resto do programa informações necessárias para a 
análise da estabilidade dos sistemas dinâmicos rotativos. Seu desenvolvimento e 
aplicabilidade são metas adicionais desta pesquisa. É importante destacar que a pa- 
rametrização da geometria possibilita também a ampliação desta biblioteca, uma vez 


que o modelo pode ser facilmente adaptado para outras configurações de mancais, 





ou até mesmo realizar adptações para outras componentes de máquinas rotativas. 


https://github.com/ross-rotordynamics/ 
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Figura 1.3: Estrutuda da biblioteca de código aberto Ross-Rotordynamics. 


Fonte: Adaptado de TIMBÓ et al.| (2019) 
1.3 Metodologia 


Inicia-se esse estudo obtendo simplificações iniciais para a equação de Navier- 
Stokes através da hipótese da utilização de um fluido laminar, incompressível e 
newtoniano, em regime permanente. Tais reduções são necessárias pois a equação 
citada possui resolução complexa, ainda que numericamente. Além disso, existem 
características geométricas que podem ser exploradas de modo a gerar simplifica- 
ções no modelo. Tradicionalmente, na literatura, Navier-Stokes é trabalhada em 


coordenadas cartesianas, negligenciando-se a curvatura apresentada na geometria 
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do problema. Com a finalidade de preservar os efeitos dessa curvatura e seguindo 


as formulações apresentadas por PINA; CARVALHO! (2006) eIANDRADE! (2008), 


as equações são trabalhadas em coordenadas cilíndricas. 


Mesmo após as considerações iniciais, ainda é necessário simplificar as equa- 
ções obtidas pois as mesmas gerariam um alto custo computacional em sua resolução. 
Em geral, os trabalhos sobre o tema realizam as simplificações seguintes com base 
em suposições razoáveis, como em (1962), que assume que a inércia e os 
termos de força de corpo nas equações de movimento são pequenos em comparação 
com os termos viscoso e de pressão. Neste trabalho, essas simplificações são rea- 


lizadas por intermédio de uma análise adimensional para o escoamento entre dois 





cilindros com uma pequena espessura de filme de óleo lubrificante. Servindo-se das 


características geométricas, especialmente para folgas muito pequenas, os termos 





das equações que descrevem o escoamento do filme fluido pudem ser ordenados por 





sua importância. Assim, para negligenciar um termo, ele deve ser quantitativamente 


Imenor que os outros. 


No limite das simplificações, é possível obter então equações mais simples 
para o campo de velocidade, sendo viável até mesmo integrá-las analiticamente, 
embora ainda dependam da pressão. As equações encontradas são agregadas devido 
à restrição da conservação de massa, modelada na continuidade, tornando-se uma 


Equação Diferencial Parcial Elíptica que tem como solução o campo de pressão. 


Por se tratar de uma equação sem solução analítica, elaborou-se um programa 
em Python a fim de resolvê-la numericamente por Diferenças Finitas Centradas. Este 
recurso, quando aplicado ao conjunto de pontos da discretização nesse tipo de pro- 
blema, resulta em um sistema linear esparso. A fim de considerar esta característica 
e reduzir o custo computacional envolvido no processo, foi aplicado o método CSC 


(Compressed Sparse Column) na matriz de coeficientes, reduzindo sua dimensão e 
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compactando as informações em suas colunas, seguido de um método adequado para 


resolver o sistema linear esparso. 


Para expressar os resultados iniciais do modelo apresentado, o campo de 
pressão é calculado para diferentes mancais, sendo feitos os devidos ajustes quanto 
a diferenças de geometria. Neste contexto, uma situação conhecida é apresentada 
e comparada aos resultados obtidos: a aproximação para um mancal curto. É 
interessante ressaltar que, na literatura, essa aproximação é feita negligenciando um 
termo do modelo já simplificado, o que possibilita realizar cálculos analíticos do 
campo de pressão. Sendo assim, este estudo visa não negligenciar esse componente 
e o campo de pressão é calculado numericamente a partir do mesmo modelo para 


mancais curtos, médios e longos, alterando apenas a geometria. 


A partir dos resultados para o campo de pressão, é possível adquirir outros 
elementos básicos para um estudo sobre mancais. Através de uma integração nu- 
mérica da pressão, utilizando a regra de Simpson composta, obtém-se as forças de 
reação do filme fluido. Tais forças, originárias da rotação do eixo, fazem com que o 
mesmo atinja uma posição de equilíbrio, calculada através de um método iterativo. 
A partir de uma posição inicial arbitrária para o centro do rotor, utiliza-se uma 
função residual que compara a carga resultante com a carga externa aplicada. O 
método iterativo apresentado faz então pequenos movimentos nesta posição, visando 


encontrar o mínimo local desta função através do método de Mínimos Quadrados. 


Por fm, supondo pequenas perturbações em torno do equilíbrio, um novo 
método é proposto a fim de calcular os coeficientes dinâmicos de rigidez e amorteci- 
mento. Normalmente, na literatura, estes coeficientes são calculados com base nas 
aproximações para o mancal curto, fundamentadas na equação de Reynolds (1.1), 
que considera a variação temporal. Dessa maneira, é possível derivar as forças em 


relação aos deslocamentos e às velocidades. Entretanto, como o presente modelo 
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analisa o caso estacionário, é apontada uma alternativa, através de uma solução por 


Mínimos Quadrados obtida pelas pseudo inversas de Moore-Penrose, que permite 





avaliar o impacto da perturbação apenas no momento em que se deseja calcular os 





coeficientes dinâmicos. Vale destacar que este é um método não usual e que torna 
a obtenção dos coeficientes computacionalmente eficaz, reduzindo o tempo gasto no 


processo. 


Intencionando fornecer ferramentas que auxiliem na elaboração de designs 
de operação, são apresentados os impactos da variação dos principais parâmetros 
de entrada na simulação de mancais cilíndricos. A parametrização escolhida para a 
geometria do problema permite também adaptar o modelo para diferentes configura- 
ções de mancais. A modelagem é aplicada a mancais elípticos (ou mancal “limão”) 
e com desgaste, analisando o impacto dessas mudanças na resposta dinâmica do 
sistema. Por fim, é realizada uma comparação entre as três geometrias, explorando 


suas alterações de comportamento para diferentes rotações. 
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2 MODELAGEM TEÓRICA 


Neste capítulo serão apresentados as características geométricas e os desen- 
volvimentos e simplificações necessários para a construção da modelagem teórica 
do presente trabalho. Este, baseia-se nas equações de Navier-Stokes e continui- 
dade. Segundo (2014), essas equações formam 
um conjunto de quatro equações diferenciais parciais não lineares que descrevem o 


escoamento. 


Por meio de uma análise adimensional e utilizando simplificações obtidas a 
partir de características geométricas, será possível obter uma equação mais simples 
de ser resolvida numericamente e que ainda assim descreve o escoamento desejado 


de forma satisfatória. 


2.1 Características geométricas e operacionais 


O escoamento do problema estudado ocorre no espaço anular entre dois ci- 


lindros de comprimento L, como ilustrado na Figura (2.1b). O cilindro externo é 





fixo, de raio R, e recebe o nome de estator, enquanto o interno, de raio R;, trata-se 
de um eixo com velocidade de rotação w, denominado rotor. É importante destacar 
que a folga radial F = R, — R; sempre será muito pequena quando comparada aos 


raios (Ro, R;) e ao comprimento (L). 





Devido à rotação, são geradas pressões no filme de óleo lubrificante e, conse- 





quentemente, forças que atuam no mesmo. Se a velocidade de rotação for constante, 





essas forças fazem com que o eixo se movimente até atingir uma determinada loca- 
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(a) Vista frontal. (b) Vista lateral. 


Figura 2.1: Geometria do problema. 


lização, chamada de posição de equilíbrio, indicada na Figura (2.1a). Nesta posição 
existe excentricidade entre os cilindros, sendo e a distância entre os centros e 5 0 


ângulo de atitude, formado entre a linha dos centros e o eixo vertical. 


Por essa razão, a descrição do rotor a partir do centro do estator varia na 
direção tangencial. Utilizando a lei dos cossenos e chamando de Rg a distância do 


centro do estator até as bordas do rotor, obtém-se a equação (2.1). 


Ro=vR?-esin?a-+ecosa, (2.1) 


3 
SR E, se S+B<0<T+B (1) 


onde a = 


-(T-0+8), e Trp<o< Ts (IT) 
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2.2 Formulação 


Inicia-se apresentando as equações de Navier-Stokes como indicado em (12.2). 


(ge tv vv)=Va. (2,2) 


onde p é a massa específica do fluido, v é o campo de velocidade cujas componentes 
são representadas por u, vew, eo = —pl+7 é o tensor de Cauchy, no qual p 


representa o campo de pressão, 7 é o tensor das tensões e I o tensor identidade. 


De acordo com |FOX:; MCDONALD; PRITCHARD| (2014), supondo a hi- 


pótese de um fluido incompressível, pode-se considerar p como uma constante in- 





dependente. Em acréscimo, o autor afirma que se o fluido escolhido for também 
newtoniano, a tensão de cisalhamento é diretamente proporcional à sua taxa de de- 


formação, sendo a viscosidade a constante de proporcionalidade, isto é, T= u(Vv). 





Além dessas hipóteses, o modelo será construído em regime permanente, descrevendo 


situações nas quais as grandezas não variam com o tempo. 


Após estas considerações iniciais, as equações de Navier-Stokes podem 
ser reescritas como em (2.3). Apesar das simplificações, uma resolução numérica 
neste momento é possível, contudo ainda com um alto custo computacional por se 
tratar de uma equação diferencial parcial não linear com uma incógnita escalar e 


outra vetorial de dimensão três. 


p(v.Vv) =-Vp+uV?v. (2.3) 


Outra importante realidade a ser considerada é a conservação de massa, des- 
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crita na equação da continuidade que, por incompressibilidade, torna-se a equação 


E3. 


a + V-(ov)=0, (2.4) 


Na literatura, uma suposição comum é que o raio de curvatura dos sólidos 
que limitam o fluido é grande em comparação à espessura do filme lubrificante, 
desprezando os efeitos da curvatura. Objetivando garantir tais efeitos e simular o 
escoamento em um espaço anular, seguindo as formulações apresentadas por 


CARVALHO! (2006) eIANDRADE! (2008), as equações serão trabalhadas em coor- 


denadas cilíndricas, como registrado a seguir: 


Direção axial, 2: 


é eu Ou e Lo du pIdu du (2.6) 
Po o ro)” ad Plrorl Or r2002 022) 


Direção radial, r: 


do du ado ué) 
á “Br a r 90 ro) 








Op 9 |19(ro)] 1Ov 20w du (2.7) 
or Or lr Or r2002 200 0272) 
Direção tangencial, 6: 
dm pu (dio tu) 
PAO Or á oz r 00 rr) 
10 9 (19uw)) 1dw 200 Sw (2.8) 
r 06 Or lr Or r2 002 1200 dz] 
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Continuidade: 


(ro) (10 Ou 
r Or ro Oz 





0. (2.9) 


2.3 Análise Adimensional 


Apesar das simplificações que reduziram as equações de Navier-Stokes, sua 
resolução ainda é complexa e de alto custo computacional. Por isso, em busca da 
compreensão da importância de cada termo da equação de Navier-Stokes, será re- 
alizada uma análise adimensional. Aproveitando as características geométricas e 
comparando a magnitude dos termos, é possível ordenar os componentes das equa- 


ções em importância e decidir quais efeitos podem ser negligenciados. 


Inicia-se considerando U e L como uma velocidade e um tamanho típicos, 





respectivamente, sendo L o comprimento com a mesma ordem de grandeza dos 
raios de cada cilindro. Como dito anteriormente, a folga radial (1) é muito pequena 
quando comparada aos raios (Ro, R;) e ao comprimento (L), sendo válida a relação 
(Ro—-R)=F<«<L. Esta característica geométrica será utilizada para classificar 
os termos da equação por ordem de importância, a exemplo do que fora citado 
anteriormente. A partir daí, pode-se desprezar alguns termos para simplificação, 
bem como identificar quais destes impactam mais a equação para futuros estudos 


Investigativos. 


O primeiro passo é separar as dimensões. As grandezas adimensionais serão 


denotadas com acento circunflexo. Já as dimensionais não serão acentuadas: 


e Velocidade: u= Ui, v=Uô, w=Uw; 


do 
e Pressão: p= Pp; 
e Variação no espaço: Az = LAZ, Ar=FAf, A(r0)= LA(?6). 


Substituindo as relações dadas acima na equação (2.6) que representa o mo- 


vimento na direção z, chega-se a sua forma adimensional em (2.10). 





“ara t lara! Le 0 


—0Pp (7 9, Lp] 1 9Uà a) 


«OU «OU Te 59) 
pIU 





ams |lt==|+ 55-25 + 555 2.10 
oz: "“lirorr |“ orr|" nm 002 “ora (2.10) 


No espaço anular a pressão apresenta maior variação na parte mais confinada 
da geometria, portanto, no parâmetro adimensional da pressão, a grandeza espacial 


, x 2) 
a ser considerada é o F e não o L, ficando da forma P = 


F? 
equação (2.10), chega-se à (2.11). 





. Reorganizando a 


pUÍ 35 + (E) 05 + a 
ER E o RT 


00) (Flda (PO da (2.11) 
“ol lpr)ror'ir/oz)) | 


UL 
Fazendo manipulações algébricas para explicitar o número de Reynolds [Re e ] 


à o. (1 
e 55 + (5 


encontra-se (2.12). Este trata-se de um conhecido adimensional que relaciona as for- 
ças de inércia do fluido (advecção) às forças de cisalhamento viscoso. Ou seja, para 
altos números de Reynolds é possível definir que o escoamento é predominantemente 


inercial, enquanto para valores baixos é predominantemente viscoso. 
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pUL 55 + (5) 05 + oo — 
u Jo lr/Co 7)” 


19 [0] (PY Ida (PN oa) 
por | Of Lp )ro Àr2)os2)) | 


F?2 


Por último, multiplicando por (=) a equação fica da forma apresentada 


em (2.13). 





R Hr 15 + (5) pOU 4 Fº w di = 
“dplêo ML) o 'lr)lro)” 

op lo | di Fe) 10% Fº à 

-5 + (r55 [195] + (77) aos + (72) 2) ea) 


2 
L? 
maneira, é possível observar em (2.13) que existem termos menores que outros. Na 





Visto que F < L, a razão (5) é pequena e menor ainda. Dessa 
teoria da lubrificação tradicional, todas as parcelas multiplicadas por essas razões 
são negligenciadas. É importante destacar que estas simplificações só serão coeren- 
tes para escoamentos viscosos, isto é, com valores baixos do número de Reynolds 
(Re). Quando este adimensional apresenta valores altos, tais parcelas podem não ser 
desprezíveis (a relação entre o número de Reynolds e tipo de escoamento na teoria 
da lubrificação é apresentada em (2010a). Ao mapear quais termos 
impactam de maneira mais significativa a equação em relação a escala de grandezas, 
observa-se, por exemplo, que o termo ds é um candidato a retornar ao modelo 
para considerar alguns efeitos inerciais. Seguindo a teoria da lubrificação tradici- 

2 
L2 
equação do movimento axial na direção axial z se reduz a (2.14). 


F 
onal, admite-se que as parcelas com os termos (5) ou são desprezíveis, a 
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op 10/00). 


De forma análoga, as equações do movimento nas direções radial r (2.7) e 


tangencial O (2.8) se transformam nas equações (2.15) e (2.16), respectivamente. 


op 
== = 21 
o [(19(7b) lop 
o? o “7000" a 


Após as simplificações, fica claro que a pressão na direção radial é constante, 


pois sua derivada é zero e a velocidade radial desapareceu. Por esse motivo, pode 





ser considerado insignificante em comparação com outras velocidades. Vale frisar 
que, havendo turbulência, v não pode ser negligenciada. Para incluir esse efeito, um 
termo candidato para retornar à equação, entre outros termos, é a velocidade radial, 


tornando-a suscetível às investigações em casos de perturbação. 


Ao retornar com as dimensões nas equações simplificadas, considerando as 
relações apresentadas anteriormente entre as grandezas e suas adimensionalizações, 


o conjunto de equações de Navier-Stokes reduze-se à (2.17) e (2.18). 


op Lo ( du. 


ceê(itm)]o am 
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2.4 Velocidades 


Com as simplificações típicas realizadas até aqui, é possível agora integrar as 
equações (2.17) e (2.18) com a finalidade de encontrar as velocidades nas direções 
z (velocidade axial u) e 6 (velocidade tangencial w). O mesmo não ocorrerá se os 


termos (5) não forem negligenciados. 


O cálculo da velocidade na direção z será realizado baseando-se na equação 
12.17). Enquanto o da velocidade na direção O será realizado a partir de (2.18). 
Reordenando-a e integrando em r, obtém-se (2.19) e (2.20). 








fa E (2 aco f (Car o CUL fr A, quo 
PES pa ao) PAR dp PR 2 
(2.20) 


Integrando novamente, chega-se as equações (2.21) e (2.22). 


1 |dpr* 
ee E +a nro | (2.21) 
“-dJj1lodp 1 CA 
v=(5 Es (inr= 5) +ar| + 6 (2020) 


onde c4, c9, c3 e cy são constantes da integração na variável r e, para encontrá-las, 


será necessário aplicar as condições de contorno. 
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2.4.1 Condições de Contorno para Velocidades 


As condições de contorno decorrem das características físicas de não desli- 
zamento e impermeabilidade adotadas. No caso analisado, apresenta-se um eixo 


interno com rotação envolto por um estator fixo. 


Devido ao princípio da aderência, sabe-se que as partículas de fluido unidas 
às superfícies sólidas adquirem as velocidades dos pontos dessas superfícies com os 


quais estão em contato. 


Por isso, como rotor e estator não se movimentam na direção axial, as con- 


dições de contorno para a componente u serão: 


Na parede do estator a velocidade tangencial também é nula. Por outro 
lado, pela rotação do rotor, na parede do mesmo a velocidade será apontada como 
o produto da velocidade angular pelo raio do rotor, acarretando as condições de 


contorno a seguir para a componente w: 
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2.4.2 Velocidades 


Aplicando as duas primeiras condições de contorno na equação da velocidade 


axial (12.21), obtém-se as constantes de integração c, e cy: 


op Rj 
Co — sea + In Ro, (2.23) 
OD: uasá A 
402 (R$ E R$) 


(2.24) 


Já as outras duas condições, aplicadas na equação da velocidade tangencial 


12.22), resultam nas constantes de integração cy e c3: 





10p 1º cR? 
e ] o pe é Za 
C4 Soo (m ; > | (2.25) 
op 


1 1 
R$- R2 


(9) 


(2.26) 


Ca — 


Substituindo as constantes provenientes da integração, descritas nas equações 
12.23) a (2.26), finalmente encontra-se as equações para as velocidades axial u (2.27) 


e tangencial w (2.28). 
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emalicr Da 
E é [a reu (E) (o) , (2.27) 
PA Ro 
1 Op 1 R$ 1 
2º 2400 Pino) tr e (mer 5) (2.28) 
wRiRo (Ro 
89 Nor): 


onde 
1 


1 1 
2 2 
+= pg [Re (no = 5) = 6 (mo 5). 


2.5 Equação da Continuidade 


Como consequência das simplificações realizadas no conjunto de equações de 
Navier-Stokes foi possível obter os perfis de velocidade do fluido. Entretanto, ambas 
as equações são descritas em função da variação da pressão. Com o propósito de 
obter esta informação, a equação da continuidade será integrada na região 


anular de interesse. 








IN Es Rm Au) dr =0. (2.29) 


Os parâmetros R, e R, do intervalo de integração podem sofrer variações 


tanto na direção tangencial quanto na axial. Desta maneira, a formulação permitirá 





incluir futuramente variações pela excentricidade e também por quaisquer caracte- 


rísticas geométricas. 
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Por propriedade, a integral da soma (2.29) poderá ser escrita como a soma 


de três integrais: 








(espe (ge pteg)a o em 


— ——— O — —  —— eee” 
I IH HI 


Para resolver as integrais de é necessário servir-se de alguns artifícios 
algébricos. Na parcela (1) é possível utilizar o teorema fundamental do cálculo 
enquanto nas parcelas (II) e (III) aplica-se a regra de Leibnitz. Detalhando cada 
termo, tem-se: 


o [o E] dr = Rov(Ro) — Rov(Ro); 


Ro Or 


Ro OW 'õ, Ro OR, OR, 
IH E — — — = ae 
db) h t3 = h id ul) do o) ag | | 


(III) /, o e, dr = E /, ru) dr — ur E do) | 

















Ponderando acerca de (1), (IN) e (II), algumas afirmações podem ser feitas: 


(1) A velocidade radial é nula em v(R,5) = O. Entretando, v(R$) £ O pois a 
origem do referencial não está no centro do rotor. 


(11) Como visto anteriormente, w(R5) = 0 e w(Rg) = wR;. Devido à excentrici- 


OR 
dade, tem-se que o Ee, 


(111) Pela condição de contorno sabe-se que u(R,) = u(Rg) = 0. 





Assim, a equação (12.30) pode ser reescrita como: 
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Ro 
E go 


Ro 
— — =), 2.31 
20 Ja, a + E» IN rudr — Rov(Ro) = 0 (281) 


A fim de calcular w(Ro) e v(R$), será considerado um ponto qualquer A 





pertencente à superfície do rotor. A posição deste ponto em relação à origem O 


do sistema de coordenadas cartesianas é chamada Ro, como representado na Figura 


Pp). 





Figura 2.2: Descrição gráfica da velocidade do ponto A. 


Fonte: Elaborada pelo autor. 


Em consequência da rotação, o ponto 4 apresenta uma velocidade v,o;, des- 
crita como a soma das componentes radial (v,ada) e tangencial (van). Sendo e, e 
eo vetores unitários do sistema de coordenadas cilíndricas, essa velocidade pode ser 


escrita como: 


Vrot = Vrad er + Van €0; (2.32) 


onde vVrad = V(Rg) € Vran = W(R). 
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Derivando o vetor posição a = Roe, em relação ao tempo, encontra-se a 


velocidade v,o: relativa a um referencial inercial. 


Urot = Ro er + Ro és. (2.33) 


Para calcular a derivada no tempo R, é necessário usar a regra da cadeia pois 
Ro é uma função de O (ângulo do ponto 4) que, por sua vez, depende do tempo t 


devido a velocidade angular. 


— 0Ro | 0R900 OR 


ec Al = a 2.34 
"oa "ooo 00” id 
Já para calcular e, utiliza-se o seguinte artifício da física cinemática: 
Cp=wxe=wlel|eg = we. (2:99) 
Substituindo (2.34) e (2.35) em (2.33), obtém-se: 
oR 
TD PR wu er +HwR,eo. (2.36) 


00 


OR 
Logo, v(R9) = wu e w(Ro) = wRo. Esses resultados podem ser substituí- 


90 
dos na integração da equação da continuidade (2.31). 
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Ro 


OR, O (Ro ORo 
— d Ro-—— = | dr — wRo-—— = 0. PARTA 
us CV da O, o 
É finalmente chega-se a: 
O (Ro O (Ro 
a d = | dr =0. 2.38 
90 mt, a 


As velocidade w e u são conhecidas e, portanto, podem ser substituídas para 


calcular as integrais, obtendo as equações (2.39) e (2.40). 


Ro PRN Es 4u 06 
1 
([Rêm R, — R$ln Ro +(R2— R)(k— 1)| — 282 (in o 4h — 5) In Ea) 


Ro 
R Rj 1 
RiR, Int) [1i+—— 2 |O 
ER o (52) ( Rr) ; 


(2.39) 
Ko 2 
/ dire = 
Ro Su Oz 
R2 o R2 = (R$ = R$) R$ o R$ RºIn (=) = fi E R$) 
a 2R$ 2) (5) “À Ro 2 
q 4H R, 


(2.40) 


Ao substituir (2.39) e (2.40) na equação da continuidade (2.38), pode-se se- 


parar os termos que dependem da pressão do lado esquerdo. Chamando de CC, a 
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parcela da equação que não depende da pressão, de €, a parte do resultado que 


multiplica o e €5 a que multiplica a encontra-se (2.41). 
096 Oz 

O op O op O 

di é g Ei E o 2.41 

o (Cr9g) “o (rg) To a 
onde 

. Ro R 1 
Co = —wR;Ro n Fa ( + na) = | ; (2.42) 


1 : di 
Ci= Fm [Rim R, — R$ln Ro e (Rj E Ro)(k E 1)| = 2R, (inn, p= 5) nm al 
(2.43) 


Ro 2 po (RR) R$ — Rá dado (MG RO) 
E qa O = Or SO so So. In (2º) Sºo “co” 
C Su r Fig 2R H R21 (=) de o (6) 2 
gHL| 5 
Ro 
(2.44) 


A equação diferencial (2.41) é uma equação diferencial parcial que descreve 





o campo de pressão p. Ainda que descrita de forma diferente devido a escolhas de 


parametrização, é equivalente à equação de Reynolds no caso permanente. Esta 





equação simplificada ainda não pode ser resolvida analiticamente, entretanto sua 





solução numérica tem um custo computacional significativamente menor em relação 
às equações de Navier-Stokes (2.3) apresentadas no início deste capítulo. Com a 
aproximação para a variação da pressão, as velocidades também poderão ser calcu- 


ladas. 


Para a resolução da equação simplificada ainda é necessário estabelecer con- 





dições de contorno na pressão que dependem das características do problema. Sem 


essas condições de contorno, a equação apresenta infinitas soluções. 
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3 SOLUÇÃO NUMÉRICA 


Através do modelo construído, foi possível transformar a equação tridimen- 
sional de Navier-Stokes com quatro incógnitas em uma equação bidimensional 
com apenas uma incógnita que descreve o campo de pressão no interior de 
componentes de máquinas rotativas. Trata-se de uma equação diferencial parcial 
de segunda ordem e, para sua solução, se faz necessário o uso de um método de 


aproximação numérica já que não possui solução analítica. 


Neste capítulo serão apresentados a discretização e o método utilizados para 
a solução numérica do modelo. A equação será resolvida pelo método de Diferen- 
ças Finitas centradas, considerando uma discretização espacial uniforme em cada 


direção, com solução aproximada para um número finito de pontos. 


3.1 Classificação da Equação 


Antes de escolher o método numérico que será utilizado na resolução da 
equação (2.41), encontrada após as simplificações realizadas no Capítulo |2| deste 
texto, é fundamental identificar a sua classificação. Segundo [FARLOW] (1993), este 
é um conceito importante porque a teoria geral e os métodos de solução geralmente se 
aplicam apenas a uma determinada classe de equações. De acordo com as definições 
apresentadas pelo autor, uma equação diferencial parcial pode ser classificada em 
seis ítens básicos: ordem, número de variáveis, linearidade, tipos de coeficientes, 


homogeneidade e tipo de equação linear. 


Prontamente, identifica-se a equação (2.41) como uma EDP de segunda or- 
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dem (ordem da derivada parcial mais alta na equação), com duas variáveis (0 e 2), 
linear e com coeficientes variáveis. Para a classificação quanto à homogeneidade 
e tipo de equação linear, é preciso analisar os coeficientes da equação trabalhada. 


FARLOW| (1993) define uma equação linear de segunda ordem em duas variáveis 


com o seguinte aspecto: 


Atz + Buy + Ciuyy + Dus + Eu, + Fu = G. (3.1) 


Comparando a equação apresentada (2.41) com o aspecto característico de 
sua classe (3.1), conclui-se que: 


ô, 
nO OC = Cs, G = 30 00 


e os demais termos são nulos (B = D= E =F=0). 


Ainda seguindo as definições apresentadas porlFARLOWI| (1993), como G 0, 


a equação também é considerada não homogênea. Por fim, a classificação 
quanto ao tipo de equação linear leva em consideração o comportamento de seus 
coeficientes. O autor indica que todas as equações que seguem o aspecto serão: 
parabólicas (A = Bº? — 4A4C = 0), hiperbólicas (A = Bº —- 44C > 0) ou elípticas 
(A = Bº -4AC < 0). Como B = 0, a classificação da equação dependerá do 
sinal dos coeficientes €; e 5 (2.44). 


Pela parametrização escolhida, estes coeficientes variam em função da geo- 
metria do problema, descrita na Seção deste trabalho, e da viscosidade (ju) do 
fluido, sendo este segundo um valor sempre positivo que não influenciará no sinal 
de C4 e Cj (2.44). Ro é o raio do estator (cilíndro externo) e Rg é a descrição 
do rotor (cilíndro interno) em função do centro do estator. Assim, pode-se consi- 
derar que R, > Rg para qualquer tamanho dos raios. Como R, e Rg, são medidas 


de comprimento em metros, sempre serão maiores do que zero. Levando em conta 
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a utilização do modelo em máquinas rotativas, os valores para os raio dificilmente 
excederão a medida de 2 metros. A fim de determinar o domínio de aplicação da 
equação e avaliar os sinais dos coeficientes, serão considerados valores dos raios no 
intervalo O < R,, Ro < 2 |m], com a folga radial igual a 1/10 de R, para garantir a 


hipótese de uma fina camada de filme fluido. 


Máximos e mínimos dos coeficientes 





—0.01 


| 
o 
[o] 
[Ni 


Coeficientes 
ó 
Era 


-— -—: C1 mínimo 

0047 — (2 mínimo 
C1 máximo 

—— C2 máximo 


—0.05 





0.00 0.25 0.50 0.75 1.00 1.25 1.50 1.75 2.00 


Ro 





Figura 3.1: Valores máximos e mínimos dos coeficientes C, (2.43) e C» (2.44) para 
medidas dos raios no intervalo O < Ro, Rg <2 ml. 


O gráfico da Figura |3.1] representa os valores máximos e mínimos que os 
coeficientes €1 e C5 podem atingir quando variados os raios no intervalo 
estabelecido. Nota-se que estes valores são sempre negativos para qualquer raio no 
domínio. Não há casos em que €4 e €5 são nulos ao mesmo tempo. Com tal 
característica, pode-se afirmar que A = Bº —- 44C < 0. Como a equação (2.41) 


descreve um fenômeno em estado estacionário e satisfaz esta propriedade, segundo 


FARLOW!| (1993), poderá ser classificada também como uma EDP elíptica. 


De acordo com ILEVEQUE| (2007), as equações elípticas fornecem problemas 


de valor de contorno em que a solução, em todos os pontos do domínio, deve ser 


oO 


determinada com base nas condições de contorno. O autor também aponta que 
isso normalmente leva a um sistema esparso muito grande de equações lineares a 
ser resolvido para os valores da incógnita em cada ponto da grade. Sendo assim, 
há a necessidade da escolha de um método numérico que considere as condições de 
contorno para a solução da equação (2.41), seguido de uma solução do sistema linear 


que otimize o processo. 


3.2 Problema de Valor de Contorno 


Uma vez determinada a classificação da equação (2.41), pode-se enunciar 
formalmente o problema a ser resolvido através de métodos numéricos. Sabe-se 
que a equação supracitada é um EDP elíptica que tem como incógnita o campo de 
pressão p, descrito em função das variáveis O e z. Na direção axial z, os valores da 
pressão são conhecidos nos extremos do cilindro, isto é, em z= 0ez=L. Já na 


direção tangencial 0, assume-se uma condição de peridiocidade visto que 0 = 0 = 27. 


Dessa maneira, simplificando a notação da equação (2.41) e fixando as con- 


dições de contorno, o problema a ser trabalhado pode ser estabelecido como: 


C1p (0, Z) 9 + Csp (0, Ap = (Co)o 
DO 0 DO b=s Bad SSD (3.2) 


p(0,2)=p(27,2) em 0<0<27 


Este é um problema de valor de contorno e, para obter uma solução, é neces- 
sário que haja uma discretização do domínio e a aplicação de um método numérico 


apropriado. 
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3.3  Discretização 


Para resolver qualquer método numérico envolvendo equações diferenciais 
parciais, o primeiro passo é a discretização da região para a qual se procura uma 
solução. Esta discretização ocorre através da definição de uma malha que, segundo 


CUMINATO; MENEGUETTE) (2013), é um conjunto finito de pontos pertencentes 


ao domínio. 








Figura 3.2: Malha retangular utilizada para discretização da equação e que repre- 
senta a planificação do espaço anular onde ocorre o escoamento do fluido. 


Fonte: Fonte: Adaptado de (ANDRADE, 2008). 


A fim de discretizar a equação a ser trabalhada, adotou-se uma malha re- 
tangular regular com N, nós na direção axial e No nós na direção tangencial, como 


indicado na Figura (3.2). Nos pontos desse domínio discretizado serão obtidas as so- 





luções referentes ao valor da pressão. O sistema de coordenadas cilíndricas adotado 





permite representar o domínio como uma malha retangular. Realizando-se um corte 
hipotético no cilindro em 6 = 0, é possível facilmente visualizar a malha resultante 
nas novas coordenadas, com a condição de que os nós de ambas as fronteiras do corte 


possuem o mesmo valor. Como constatado através da análise adimensional feita na 
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Seção deste texto, a pressão na direção radial é constante. Por esta razão, não 


será realizada uma discretização em r. 


Dada a malha regular, o espaçamento entre os N, pontos é nomeado Az, 
enquanto o espaçamento na direção tangencial é chamado A6, e são definidos da 


seguinte forma: 


L 27 
AO = | 
N-1 No-1 








(3.3) 


3.4 Método das Diferenças Finitas Centradas 


Com o objetivo de encontrar a solução numérica do problema, será utilizado o 
método das diferenças finitas centradas. explica que a 
ideia básica de difereças finitas é transformar a resolução de uma equação diferencial 
em um sistema de equações algébricas no qual as derivadas são substituídas por 


diferenças que usam apenas os valores numéricos da função. 


Pelo método das diferenças finitas centradas, a pressão em cada nó se relaci- 
ona com seus quatro nós vizinhos. Dessa maneira, o sistema de equações algébricas 


terá dimensão dada pelo produto de N, por No. À molécula computacional apresen- 





tada na Figura (3.3) mostra a tradução gráfica da relação de dependência existente 


entre o valor da pressão no ponto (%,)) e seus vizinhos. 


No método numérico adotado, as derivadas relativas à pressão são aproxima- 
das pela diferença dos valores nos pontos da malha, dividida pela distância entre 


eles. 


do 


(+1,)) 





--—-—. 


(1,]-1) 


(1-1, )) 
Ae 


Figura 3.3: Tradução gráfica do método das diferenças finitas centradas. 


Fonte: Adaptado de (ANDRADE/ 2008). 





Ê 
i 
! 
+——————————— 


op 
oz 


- aproxima-se as derivadas 


Assim, fazendo F, = Ci(z, ya e Fo = C5(2,6) 


por: 


0F, FE-FW FP, EN-FS 
a e EE 





— = — a 
096 2A0 Oz Zina E) 
Também o termo independente da pressão: 
Co ee. 
0 240 * No 


Escrevendo as parcelas da equação diferencial da forma discretizada, obtém- 


Se. 


4 


0) Op, 1 Pijj41 — did) (Pia = a] 

ô, Op| 1 Pi+1,j — did) (2 = Pts )] 

ER [Ed = ne [er ( As Coçi-1,)) Az , (3.6) 
Om — Cos) — Co(i-1) 

39 (00) — AD 


O esquema de discretização (3.6) usado na EDP elíptica (2.41) gera uma 


equação algébrica para cada nó (1,9), indicada na equação (3.7). 


(Co(i-1,5)) (Cr(ij-1)) 
Pis ADO + Pis ADO = 
E upsCina (O quromes 
ho ç TO) é ge ne (1.9) E a) 4 (3.7) 


(C1(i,9)) (Cox) 1 
Piij+1 nr É Pix, e] > AO Cowt(id) — Cow(ij-1) 








3.4.1 Condições de Contorno para Pressão 


Como visto na Seção os valores da pressão nos extremos do cilindro são 





conhecidos e representam as condições de contorno na direção axial. Estas condições 


serão aplicadas na equação discretizada (3.7). 


p(z=0)=p(1,5) = Pim, 
plz = L) ED ad) E E ui 


dO 


Na direção tangencial, as fronteiras da direita e da esquerda são apenas um 
artifício de corte para facilitar a implementação numérica. Como não representam 


contorno físico, esses pontos assumem uma condição de periodicidade. 


Vale lembrar que para as aproximações por diferenças finitas centradas, nos 


nós localizados em 6 = 27 considera-se que: 


op Piz 7 PinNç- 
00 240 di 


IN 


3.4.2 Resolução do Sistema de Equações 


A equação (8.7), quando aplicada ao conjunto de pontos da malha, resulta 
em um sistema linear com (N, - N9) equações. Tal sistema linear pode ser escrito 


da forma matricial (3.9): 


MI ol = fl, (3.9) 


onde a matriz M é composta pelos coeficientes que multiplicam as pressões em cada 


ponto do domínio discreto, p é o vetor de incógnitas e f representa o lado direito 





da equação, isto é, os termos independentes da pressão. 


A matriz M obtida possui uma estrutura de banda a menos da condição de 


periodicidade, como representado na Figura (83.4). Isso significa que, dependendo 
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Representação gráfica da matriz M 
para N, - 8e N,- 32 


Ô 50 100 150 200 250 


Figura 3.4: Representação esquemática da matriz M. 


Fonte: Elaborada pelo autor. 


do método escolhido, o problema pode se tornar menos esparso, exigindo maior 


memória computacional. 


Inicialmente, utilizou-se o método linalg.solve da biblioteca NumPy, por 
(2006), para resolver o sistema linear (3.9). Este método é capaz de 
resolver uma equação de matriz linear ou um sistema de equações escalares lineares, 
calculando a solução da equação da matriz linear quando bem determinada, como 
é o caso. Entretanto, a fm de considerar a característica de esparsidade da matriz 
M, esse método foi substituido por outros dois da biblioteca SciPy, por [VIRTANEN] 
(2020): primeiro o sparse.csc matrix, que aplica o método CSC (Compressed 
Sparse Column) na matriz M, reduzindo sua dimensão e compactando as informa- 
ções em suas colunas; seguido do segundo método, sparse.linalg.spsolve, que resolve 
o sistema linear esparso. A alteração na estratégia utilizada para a resolução do 
sistema linear possibilitou uma significativa redução do tempo computacional envol- 


vido no processo. Para uma discretização de N, = 16 e Ng, = 264, a mesma solução 


o 


do sistema linear foi obtida em 0.8125 segundos, enquanto levava 4.78125 segundos 
com o primeiro método. Ambos os resultados foram obtidos através de simulações 
realizadas em um notebook com processador Intel Core ài7 8565U O 1.80GHz, com 


memória RAM de 8GB Q 2400/2660 MHz e disco HD de 1 TB Q 5400 rpm. 


3.5 Convergência de Malha 


Objetivando testar a convergência do método e escolher o melhor tamanho 
para a malha a ser trabalhada, foram gerados resultados para um número crescente 


de pontos nas direções axial (N,) e tangencial (Ng). 


O teste foi composto por duas etapas: 


e Primeira etapa: fixou-se o número de pontos na direção tangencial e a con- 


vergência da malha foi verificada a partir do aumento gradativo de N,. 


e Segunda etapa: encontrado o melhor resultado para o número de pontos 
na direção axial, este valor é fixado e a convergência da malha é verificada a 


partir da variação de No. 


A verificação da convergência ocorreu através do cálculo da matriz de pressão, 
observando-se o lugar geométrico onde a pressão apresenta o maior valor. Desta 
maneira foi possível comparar as duas respostas da variação da pressão para os 


mesmos pontos em cada etapa. 


À medida que o campo de pressão era calculado para um determinado ta- 


manho da malha, também era calculado para um um número de pontos maior, 


o8 


verificando o erro relativo gerado e tendo como referência os valores obtidos pela 


maior malha. 


O teste foi realizado para a uma geometria de dois cilindros excêntricos, cujos 


parâmetros de entrada encontram-se na Tabela (83.1). 


Elemento Valor 
0.2 [m] 

0.2002 [m| 
0.0001 [ml 





Raio do rotor (R;) 

Raio do estator (R$) 
Excentricidade (e) 
Comprimento (L) | 0.4002 [m| 


0.015 [Pa.s] 
10.472 [rad/'s| 


Viscosidade (ju 
Rotação (w 


Tabela 3.1: Parâmetros utilizados para a análise de convergência. 


Os resultados da primeira etapa são apresentados no gráfico da Figura (3.5). 


O número de pontos na direção tangencial foi fixado em N, = 132 enquanto o número 





de pontos na direção axial variou em 4 < N, < 80. Percebe-se que o erro relativo a 
partir de N, = 16 torna-se menor que 0.1%, sendo portanto, o valor escolhido para 


a discretização na diração Z. 


Por sua vez, os resultados da segunda etapa encontram-se na Figura (13.6) e 
descrevem os erros relativos obtidos a partir da análise de convergência para Np. Já 
escolhido o valor para N,, a discretização na direção tangencial variou em 4 < Ng < 
200. Foi escolhido No = 132 pois, a partir desse valor, o erro relativo apresenta 


maior estabilidade e valor na ordem de 1072. 


Para a discretização escolhida, com N, = 132, N, = 16 e utilizando os 
mesmos dados da Tabela (3.1), obtém-se o comportamento da pressão indicado na 


Figura (3.7). 
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Variação de NZ para o cálculo da pressão 


0.08 


0.06 


0.04 


Erro relativo 


0.02 


0.00 
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Número de pontos NZ 


Figura 3.5: Convergência da malha para N,. 


3.5.1 Malha Quadrada 


A metodologia estabelecida para avaliar a convergência de malha gera resul- 
tados promissores. Entretanto, ao fixar uma quantidade de pontos em determinada 


direção e aumentar gradativamente a discretização na outra, há uma intensa dis- 





torção da malha, fazendo A6 cada vez mais diferente de Az. Para a malha já 


convergida, Az = 0.027 |m| enquanto Ad = 0.048 |rad)]. 


Dessa maneira, fixando A6 = Az, é possível descrever N, em função de N,. 
Com a discretização N, = 16, Aô será aproximadamente igual a Az para Ng = 236. 
Considerando essa nova malha quadrada e comparando com a malha convergida 
encontrada pela metodologia anterior, no lugar geométrico onde a pressão apresenta 
o maior valor, obtém-se um erro relativo na ordem de 107%, tendo como referência a 
malha quadrada. Deste resultado conclui-se que a solução numérica, além de conver- 
gir para o aumento da malha, também se comporta de maneira satisfatória mesmo 


que o conjunto de pontos da discretização apresente distorções em seu espaçamento. 
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Variação de N8 para o cálculo da pressão 


0.40 


0.35 


Erro relativo 
o o oo 
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[6 ,| [ao] 6 ,| [eo] 


o 
o 
[es] 


0 25 50 715 100 125 150 L75 200 
Número de pontos Ng 


Figura 3.6: Convergência da malha para Na. 
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Figura 3.7: Campo de pressão obtido para a geometria de dois cilindros excêntricos. 
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4 MANCAIS HIDRODINÂMICOS 


A escolha da parametrização do modelo apresentado neste trabalho permite 
fácil adaptação a diferentes geometrias e condições operacionais de componentes 
de máquinas rotativas. Uma vez obtida a solução numérica do modelo após as 
simplificações realizadas ao longo deste texto, é possível aplicá-lo a contextos mais 
específicos, permitindo uma maior exploração dos resultados. Por esta razão, este 
capítulo apresentará conceitos e mudanças necessárias a ffm de utilizar a modelagem 


proposta para mancais hidrodinâmicos. 


Um mancal pode ser definido como uma estrutura mecânica fixa que funciona 
como suporte de apoio para eixos e rolamentos de máquinas rotativas. Em geral 
é feito de ferro fundido ou aço e classificado de acordo com o tipo de operação, 
movimentos permitidos ou direções das cargas aplicadas nas peças rotativas. Seu 
principal objetivo é sustentar o eixo, restringindo o movimento à posição desejada 
e diminuindo o atrito entre as partes móveis. Existem vários tipos de mancais, 
entretanto, neste texto serão considerados apenas os mancais hidrodinâmicos, isto 


é, mancais que possuem uma camada de filme lubrificante entre as estruturas rígidas. 


De acordo com MEGGIOLARO! (2008), os mancais hidrodinâmicos são nor- 


malmente utilizados em máquinas rotativas de grande porte por possuírem capaci- 
dade de carga satisfatória para esta condição. Além disso, características de rigidez e 
amortecimento ajudam na condição de estabilidade dinâmica dos sistemas rotativos, 
o que justifica a sua grande usabilidade, como em motores de combustão interna, 
turbocompressores, turbogeradores, e máquinas rotativas em geral. A partir do 
movimento de rotação do eixo são geradas pressões no lubrificante que separam as 


superfícies do mancal e do eixo por uma fina camada de filme fluido em formato 
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de cunha. Previne-se assim o contato entre as superfícies metálicas, prolongando a 


vida útil dos equipamentos mecânicos. 





Figura 4.1: Desenho esquemático de um mancal. 


Fonte: Elaborada pelo autor. 


Para o caso de mancais hidrodinâmicos de estrutura simples, o escoamento 
também ocorre no espaço anular entre dois cilindros. A descrição da geometria 
nesta situação é muito semelhante à modelagem apresentada na Seção deste 
texto. Com o propósito de compatibilizar com a literatura, o ângulo de trabalho 6 


é contado a partir do maior espaçamento, como indica a Figura 4.1 


4.1 Condição de Cavitação 


Através das soluções numéricas apresentadas no Capítulo |3] deste texto, é 
possível perceber que a distribuição de pressão está relacionada à variação na es- 


pessura do filme fluido. A cunha convergente para o menor espaçamento incide 





num campo de pressão positivo, enquanto a cunha divergente reflete um campo de 


pressão negativo, como representado nas Figuras (a-b). RAMOS) (2019) indica 
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que a existência dessa cunha divergente de óleo no mancal hidrodinâmico torna a 
cavitação um fenômeno característico neste tipo de componete. Cavitação, segundo 
(1979), pode ser definida como o fenômeno que descreve a 
descontinuidade de um fluido pela existência de gases ou vapor. Os autores afirmam 
que a maioria dos fluidos, apesar de suportarem elevadas pressões positivas, não são 


capazes de suportar elevadas pressões negativas, provocando sua ruptura. 








ij . “ e “ o. - e Lubrificante 





Direção de movimento da superfície 


(a) Cunha do filme de óleo lubrificante. 
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(b) Distribuição de pressão sem considerar o fenômeno da cavitação. 


Figura 4.2: Fenômeno da cavitação. 


Fonte: Adaptado de (2019) 
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De acordo com RAMOS)! (2019), existem duas formas de cavitação nos man- 


cais hidrodinâmicos: a cavitação por gás e a cavitação por vapor. À primeira 


CDx 





caracterizada por gases do ambiente externo que migram para o fluido e, devido à 
queda de pressão, atingem o seu limite de saturação, formando fases gasosas. Na 
segunda forma, o vapor é gerado dentro do próprio fluido lubrificante quando a 


pressão é inferior a sua pressão de vapor, formando também a fase gasosa. 


Como consequência da cavitação, o autor apresenta o fato da distribuição de 





pressão manter-se positiva na parte convergente da cunha de óleo, entretanto, na 


região divergente a pressão torna-se constante e igual à pressão de vapor do fluido. 


É importante destacar que esta mudança de comportamento da pressão em 
decorrência da cavitação não se inicia necessariamente no ponto de menor espessura 


do espaço anular. Diversos estudos buscam estabelecer as condições de contorno ade- 





quadas para descrever o início da cavitação no fluido, dentre eles segundo ISHIDA: 


AMAMOTO! (2012), destacam-se: 


e Condição de Sommerfeld: Não considera a existência de cavitação. Na 
posição de equilíbrio, a distribuição de pressão é antissimétrica em relação ao 
ponto O = 7, sendo positiva na região de0 < 6 < 7 e negativa em 7 < 0 < 27. 


Essa distribuição de pressão é mantida apenas em casos onde a pressão é muito 


pequena (Figura [4.3a). 


e Condição de Gumbel: Utilizando o argumento de que podem ocorrer a 
evaporação do lubrificante e o fluxo de ar axial de ambas as extremidades, 


considera-se a pressão na região 7 < 6 < 27 como quase zero (isto é, a pressão 





atmosférica). Define-se então p = O em toda a região divergente (Figura |4.3b). 


e Condição de Reynolds: Esta condição considera que após entrar na região 


de cavitação, a pressão se torna constante e igual à pressão de cavitação do 
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fluido, ou seja, o gradiente de pressão é nulo nesta área (Figura |4.3c). 





(a) Condição de Sommerfeld. (lb) Condição de Gumbel. — (c) Condição de Reynolds. 


Figura 4.3: Condições de contorno para a cavitação. 


Fonte: Adaptado de |ISHIDA; YAMAMOTO| (2012) 


ISHIDA; YAMAMOTO|(2012) afirmam que a condição de Reynolds expressa 


a situação prática com mais precisão, mas a condição de Gumbel é amplamente usada 


por causa de sua simplicidade. Além disso, de acordo com |ISAN'TOS! (1995), embora 


viole a conservação de massa, esta última condição apresenta erros aceitáveis nos 





parâmetros globais do mancal. Por estas razões, o presente texto adotará a condição 


de contorno de Gumbel para descrever o fenômeno da cavitação. 


4.2 Aproximação para um Mancal Curto 


Na bibliografia sobre mancais hidrodinâmicos, vários trabalhos utilizam a 
equação de Reynolds (1.1), após uma série de simplificações, para encontrar o com- 
portamento da pressão em mancais. Entretanto, como trata-se de uma equação que 
não possui solução analítica, usam o artifício de aproximar a equação para os casos 
de mancais curtos (L/D > 0) e infinitamente longos (L/D — c0). Deste modo, 


negligencia-se uma das parcelas da equação, sendo possível encontrar modelos redu- 
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zidos que podem ser resolvidos analiticamente. 


De acordo com SAN ANDRES, (2010d), a maioria dos mancais modernos 


em aplicações de turbomáquinas de alto desempenho tem a relação L/D pequena, 





raramente excedendo a unidade. Portanto, o modelo para um mancal infinitamente 
longo é de interesse atual limitado. Em contrapartida, o autor indica que o modelo 
curto fornece resultados precisos para mancais cilíndricos com a razão L/D < 0.5, 
sendo amplamente usado para estimativas rápidas das características de desempenho 


das forças estática e dinâmica do mancal. 





Com a finalidade de verificar as simplificações realizadas neste trabalho e 
comparar as adaptações da geometria para o contexto supracitado neste capítulo, 
foram realizadas comparações entre a aproximação para um mancal curto e o pre- 


sente modelo. 


Neste contexto, o comprimento do mancal é considerado muito pequeno e, 
segundo |[ISHIDA: YAMAMOTO| (2012), a variação de pressão na direção z pode 
ser considerada muito maior que na direção x, ou seja, Op/0x < Op/0z. Assim, o 


primeiro termo da equação de Reynolds (1.1) é negligenciado, reduzindo-se a: 


o [500 9h 
az [h 3 = Gu (04 Ui) dg + 2 (4.1) 


A equação (4.1) possui solução analítica. Fazendo os devidos ajustes para o sistema 
de coordenadas adotado neste trabalho, obtém-se então uma fórmula que descreve 
o comportamento da pressão no mancal curto, como indicado em (4.2). Os autores 
apresentam a equação com as variáveis temporais, mas indicam que, para o caso 
estático, estas devem ser igualadas a zero. 


—3uew sin O | a) (4.2) 
curto — Do ANDI AB o “ C& ls é 
Porto (py RJ? (1 + ecoso) 2) 4 


€ / LA . e 
onde ce=> ——— é a razão de excentricidade. 


io 
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4.2.1 Comparação 


Devido à grande utilização do mancal curto no cenário atual, comparar o 
modelo proposto neste trabalho com a aproximação utilizada pela literatura para 
este tipo de mancal significa verificar a concordância da solução numérica com grande 
parte da bibliografia sobre o tema. Para esta finalidade será considerada a simulação 
numérica cujos dados estão apresentados na Tabela sendo explorada a relação 


ED: 


Elemento Valor 
0.2 [m|] 
0.0001 [ml 
0.25: F ml] 





Excentricidade (e 


0.015 [Pa.s] 
10.472 [rad/s| 


Viscosidade (ju 
Rotação (w 


| 
Ângulo de atitude (8) | 7/12 [rad] 


Tabela 4.1: Parâmetros utilizados para a comparação com a aproximação do mancal 
curto. 


A classificação do mancal quanto ao seu comprimento, é apresentada por 
(1990). Segundo o autor, um mancal pode ser considerado curto se 
a relação L/D for menor ou igual a 1/8 = 0.125. Por outro lado, 
(2010d) informa que são obtidos valores aceitáveis com a razão L/D < 0.5. Ou- 
tros autores distintos ainda discordam sobre a definição deste parâmetro. Desse 
modo, torna-se importante analisar até que ponto a aproximação da Equação (4.2) 


é aceitável. 


Ip o Dsmrtálos 


Erro relativo = 
| einio | | o 


(4.3) 
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O gráfico da Figura[4.4a]apresenta o erro relativo (4.3) da matriz numérica de 


pressão, tendo como referência a solução aproximada (4.2). Variando o parâmetro 








L/D no intervalo entre 0 e 1 é possível perceber que, como esperado, o erro cresce a 
medida que o mancal torna-se maior. Isso acontece devido à parcela negligenciada 
na equação de Reynolds Para valores de L/D < 0.5 o erro relativo aproxima-se 


de 0.1, validando a afirmação feita poriSAN ANDRES! (2010d). Entretanto, os erros 


atingem um valor menor que 0.02 a partir de L/D < 0.2, como indica a Figura [4.4] 


Dando destaque ao intervalo O < L/D < 0.2, o gráfico da Figura [4.4b) expli- 
cita que não é monótona a relação de que quanto menor o comprimento do man- 
cal, melhor será a aproximação. Também é possível observar o comportamento 
da pressão no meio do mancal para diferentes comprimentos na Figura |4.5, Com 


L/D = 0.112 obtém-se o menor valor do erro relativo (menor que 0.009). Porém, 





a esquerda desta razão, percebe-se um novo aumento do erro, indicando que existe 
um intervalo no qual a aproximação é mais adequada dentro de uma determinada 
tolerância. Estabelecendo-se o critério de erros menores ou iguais a 0.01, o intervalo 


no qual a aproximação para o mancal curto é indicada será 0.088 < L/D < 0.134. 


Adotando a proporção L/D = 0.125 indicada por FRÊNE et al.| (1990), o 


modelo proposto neste trabalho apresenta um erro relativo E = 0.009 quando com- 
parado a solução para mancais curtos de (2012). O gráfico 
comparativo é ilustrado na Figura |4.0] Este resultado preliminar é extremamente 
relevante uma vez que esta aproximação é utilizada em diversas simulações para 
mancais, indicando que o modelo concorda com a literatura. Além disso, as aná- 
lises do comprimento do mancal reforçam o fato de que a completude da solução 
numérica ganha importância a medida que o mancal torna-se maior, dado que as 


simplificações realizadas na literatura não se comportam bem nestas situações. 
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Comparação do erro da aproximação da pressão para diferentes razões L/D 


0.25 


Erros 
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LD 


(a) Erro relativo da pressão no mancal curto para O < L/D < 1. 


Comparação do erro da aproximação da pressão para diferentes razões L/D 
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Figura 4.4: Erros da pressão para diferentes valores da razão L/D. 
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Figura 4.5: Pressão ao longo de O no meio do mancal para diferentes valores de 
L/D. Os valores pontilhados representam a solução numérica do presente trabalho, 
enquanto as linhas contínuas indicam os resultados da aproximação apresentada por 


ISHIDA; YAMAMOTO| (BOT3). 


4.3 Forças do Filme Fluido 


Uma vez obtido o campo de pressão, agora é possível obter as forças de reação 
que atuam no eixo. Este estudo é importante pois tais forças podem levar o rotor 
à instabilidade, sendo fundamental para análises de sistemas dinâmicos. De acordo 
com (2010d), a integração do campo de pressão na superfície do 
mancal gera os valores das forças do filme de óleo nas direções radial N (oposta à 


excentricidade) e tangencial T indicadas na Figura [4.1] dadas por: 





N L per Icos6 
r|= / / » [esa] dodz. (4.4) 


Como o comportamento da pressão é obtido numericamente, as integrais (4.4) 


(1 
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Figura 4.6: Pressão ao longo de 0 no meio do mancal curto com L/D = 0.125. 


também necessitam de um método numérico para serem resolvidas. Para isto foi 


escolhida a regra de Simpson composta aplicada através do método integrate.simps 


da biblioteca SciPy, por [VIRTANEN et al. (2020). 


Para o caso da aproximação para um mancal curto, substituindo o compor- 


tamento da pressão dado em (4.2), é possível realizar analiticamente as integrais 


e obter uma equação para as forças, dadas por ISHIDA:; YAMAMOTO! (2012) no 


contexto estacionário como: 


1 RR. Sºb | dew 
N=u(r>) —|nal 
Z la = dt; Fr (1— e?) 
1 R; a Ee TEW 
| 
2 sd dp ifh=o e 


As forças resolvidas numericamente quando comparadas às forças (4.5) do 


(4.5) 
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mancal curto, na mesma simulação numérica apresentada na Figura |4.0l apresentam 
erros relativos, tendo como referência a aproximação para o mancal curto, de Ey = 


0.02 para a força radial e Er = 0.005 para a força tangencial. 


Através da transformação (4.6) das componetes N e T também é possível 


obter as forças f; e fy no sistema de coordenadas cartesianas: 


fu =Tcos(5) — Nsin (5), 
fy = Tsin(B) + Ncos(5). (4.6) 


4.3.1 Posição de Equilíbrio 


Pela definição de mancais hidrodinâmicos apresentada no início deste capí- 
tulo, sabe-se que estes são estruturas desenvolvidas para sustentar um eixo com 
rotação. Este eixo possui uma carga de peso vertical W como indicado na Figura 


De acordo com |FRISWELL et al.| (2010), é muito comum que exista uma carga 


resultante, também vertical, devido a este peso do rotor. 


Como indicado por SAN ANDRES! (2010d), em condição de equilíbrio, isto 


é, quando o rotor atinge o deslocamento excêntrico e com ângulo de atitude 5, a 
força resultante vertical se equilibra com o carregamento externo aplicado W na 


velocidade de rotação w. Assim, as equações de equilíbrio estático serão: 


Tn =W. (4.7) 


Conhecendo o carregamento externo W, é possível através de um método iterativo, 
chegar à posicação de equilíbrio. A partir de uma posição inicial para o centro do 
rotor, utiliza-se uma função residual que tem como imagem a componente f, das 


forças e a diferença entre a componente f, e a carga W. O método iterativo então 


(3 


faz pequenos movimentos no centro do rotor, visando encontrar o mínimo local desta 
função. Para tanto, foi empregada a ferramenta optimize.least squares, também da 
biblioteca SciPy de [VIRTANEN et al (2020). Para fins de otimização, as variações 
na posição do centro do rotor são limitadas ao quarto quadrante para uma rotação 


do sentido anti horário. 





Figura 4.7: Esquema do método iterativo para encontrar a posição de equilíbrio no 
mancal. 


Fonte: Elaborada pelo autor. 


Utilizando o método apresentado na simulação descrita na Tabela [4.2] encontra- 


se, satisfatoriamente, fr = |f,— W| = 10º[N]. 


Elemento Valor 
0.2 [m] 
0.0001 [m| 
0.125 





Raio do rotor (R;) 
Folga (HF) 
(L/D) 

Carga (W) | 50 [N] 


0.015 [Pa.s] 


Viscosidade (ju 
( 10.472 [rad/s| 


Rotação (w 


Tabela 4.2: Simulação para um mancal curto sendo fornecida a carga externa apli- 
cada. 
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De acordo com [SAN ANDRES! (2010d), as características de desempenho 


estático do mancal são relacionadas em um adimensional conhecido como Número 


de Sommerfeld (5): 


s- uwLR? 


—— 4. 
nW F?2 ) ( 8) 


onde EF = R, — R, é a folga radial. 


Este adimensional é diretamente ligado à posição de equilíbrio. Segundo 
(2012), mesmo que a rotação e a viscosidade mudem, a 
posição de equilíbrio não muda se o número de Sommerfeld for o mesmo. A Figura 
[É.8l ilustra essa relação. Quando não há rotação (ou seja, S = 0), o rotor localiza-se 


na parte inferior do mancal. A medida que a velocidade de rotação aumenta (e 





consequentemente S também torna-se maior), o centro do rotor flutua para cima, 
percorrendo um caminho em formato de semicírculo até chegar ao centro do estator, 


quando a velocidade de rotação se torna extremamente alta (ou seja, S = 00). 





S=0,€=1 


Figura 4.8: Caminho percorrido pelo centro do rotor relacionado ao Número de 
Sommerfeld. 


Fonte: Adaptado de!FRISWELL et al. (2010). 


To 


4.4 Coeficientes Dinâmicos 


As forças do filme de fluido podem ser consideradas funções gerais dos des- 
locamentos e velocidades do centro do rotor. Supondo pequenas perturbações em 


torno da posição de equilíbrio, conforme SAN ANDRES, (2010e), permite-se então 


expressar tais forças como uma expansão em Série de Taylor: 

















AN A 
fe ao É + SE ny + Ss + Eng + Gs 
oh fu ng Pons Tung Lungs a; 
fy = Jo É Sr Ag ad Ay VA Sã a Ay V+ 5% Rd (4.9) 


A partir desta expansão, o autor também define os coeficientes de rigidez 
(k;;), amortecimento (c;;) e força de inércia (m;;) como: 
o fi o fi o fi 


p= ss Cij = — 7) Mij = — 5) 1] EU. 4.10 
j 0; j 05 j 0) J y ( 











De acordo com (2017), os coeficientes dinâmicos de um mancal são 
componentes necessários na análise da estabilidade linear e resposta de sistemas di- 
nâmicos rotativos. Segundo o autor, a estabilidade do sistema é uma consideração 
importante para a maioria das máquinas rotativas de alta velocidade, e os coefici- 
entes dinâmicos de um mancal são de grande importância na dinâmica e na análise 


dessa estabilidade. 


Em (4.10), os elementos nos quais 1 = 3 são chamados de termos diretos, 
enquanto para à * j são denominados como termos cruzados. Os coeficientes de 
inércia do fluido (mi;), segundo [SAN ANDRES) (20106), são de grande importância 
em mancais de escoamento superlaminar e turbulento, principalmente para líquidos 
densos. No entanto, seu efeito na resposta dinâmica dos sistemas de mancais é 
importante apenas em altas frequências de excitação. Por esta razão, este texto não 


considerará estes coeficientes. 
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Segundo |[LUND) (1987), a característica carga-deslocamento de um mancal 


é evidentemente não linear, sendo a definição dos coeficientes de rigidez e amor- 





tecimento apresentada uma solução linearizada. O autor afirma que essa solução 





concorda bem para o caso de perturbações de pequena amplitude, mas é incapaz de 


prever as órbitas concorrentes com amplitudes maiores. 


4.4.1 Perturbação no centro do rotor 


Todas a análises realizadas até agora consideram o rotor na posição de equi- 
líbrio estático. Entretanto, para obter os valores dos coeficientes dinâmicos é neces- 
sário analisar o comportamento das forças por meio de uma pequena perturbação 


em torno deste equilíbrio. 





Considere novamente um ponto 4 qualquer pertencente à superfície do rotor. 





Quando uma velocidade de perturbação é imposta, ou seja, o centro do rotor é 
movido em uma pequena amplitude, 4 também se desloca, adquirindo uma nova 


velocidade, descrita em (4.11) e exemplificada na Figura [L.9) 
Utot = Vrot + Up; (4.11) 


onde vo é a velocidade do ponto devido à rotação (2.33), descrita na Seção 


deste texto, e v, é a velocidade adquirida pela perturbação. 


Deslocando o centro do rotor horizontalmente, sua nova abscissa será dada 
por x = x0 + xp sin(wpt), enquanto ao deslocar verticalmente, a nova ordenada será 
y = Yo + ysin(wst). Os termos com O subscrito indicam os valores na posição 
de equilíbrio, enquanto aqueles com a letra p estão relacionados à perturbação, 
wp é a frequência de excitação e t é o tempo. À transformação entre o sistema de 


coordenadas cartesiano e o cilíndrico é (considerando que e, e eg são vetores unitários 


vt 





Figura 4.9: Descrição gráfica da velocidade do ponto 4 com uma perturbação na 
direção x. 


Fonte: Elaborada pelo autor. 


do sistema de coordenadas cilíndricas): 


er = ercos(0) — eg sin(6), 


ey = er sin(0) + eg cos(0). (4.12) 


Assim, para uma perturbação em cada direção, v, é dada por: 


lar cos ul o 


(vp)y = WpYp COS (wpt ey. (4.13) 
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É então, 


tp E + wWpZXp COS (wpt) COS o) er + (wRg — wpxp cos (wpt) sin (0))eo, 
E —————— 


“08 
(Ro) Md 
UVlito)x 
ORo 
Ut) DT + wpyp COS (wpt) sin (0) [er + (wRo — wpyp COS (Wnt) cos (0) ) eo. 
a 
PD qt w(Ro)y 


(4.14) 


Substituindo as velocidades radiais (v(Ro)r e v(Rg),) e tangenciais (w(Ro)s 
e w(Ro)y) para as perturbações nas direções x e y (4.14) na análise da equação 


de continuidade (2.31), encontra-se, para ambos os casos, uma nova parcela deri- 





vada em relação a 0, independentemente da pressão. Isso significa que, na equação 
simplificada (2.41), todo o lado esquerdo da equação será mantido. Entretanto é 


necessário adicionar uma nova parcela ao coeficiente Co: 


(Co): = Co + [Rowpxp cos (wrt) sin (6)], 
(Co)y = Co + [Rowpyp cos (wpt) cos (6)]. (4.15) 


4.4.2 Coeficientes de Rigidez e Amortecimento 


Ao incluir os efeitos da perturbação na equação simplificada (2.41), torna-se 
possível reescrever a expansão da Série de Taylor das forças (4.9) em função do 


instante de tempo t para cada direção da perturbação. 


(Set), = Loo + her A xft] + cor Stft), 
ul e Rg || Ga] 

(fot | Aúft] 
Cl 


fylt y Juo + hyyylt] + coyAgft] 


y = FER E AN t| E Cry 


) 

p 4 
) 
) | 


(4.16) 
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Escrevendo (4.16) na forma matricial: 


e Na direção 1: 





AcO)| AilO| (Le [0Do — Fo (y [OD — fo 
Ar) as | ba tel | De fo Do fa 
LAc|N 1) AN) (EIN = De — fes (IN = We — fool 
o º (4.17) 
e Na direção y: 
Aylo] Ay(O] (Sa 10) dm (1, [0], — Ju 
Ay] Ay] ps Ed o (o Dy — fxo Cy Dy — fyo 
AyN 1] AjN-1 GEN = My — fes (Fo = Dy — fa 
o º (4.18) 
Ou, de forma a simplificar a notação: 
afertel=r e afetos qu 


A fim de calcular os valores dos coeficientes será utilizado o método de Míni- 





mos Quadrados. Para tanto, utilizam-se as pseudo inversas de Moore-Penrose das 
matrizes A, e A,. 


AL=(ATA,) AT ec AS=(ATA) AT 


uy: 


(4.20) 
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É, finalmente, chega-se aos coeficientes de rigidez e amortecimento: 


Ky Cy 
dá dá | K xys Cy 
estator 7 o? 
/ ) 


, Rd 
K yX ' Cy a Pd 


Figura 4.10: Coeficientes de rigidez e amortecimento do mancal. 


Fonte: Adaptada de (BOMPOS D. A.; NIKOLAKOPOULOS 2014) 








É importante destacar que, em geral, os trabalhos relacionados não conside- 
ram a equação de Reynolds no caso estacionário, mas sim com a variação no 
tempo. Ao realizar as aproximações para o mancal curto, é possível então derivar as 
forças em relação aos deslocamentos (Ax e Ay) e às velocidades (Ax e Ay). Con- 
tudo, a alternativa encontrada neste trabalho através de mínimos quadrados permite 


continuar considerando o modelo no caso estacionário e avaliar o impacto da per- 





turbação apenas no momento em que se deseja calcular os coeficientes dinâmicos. 





Este fato faz com que a obtenção dos coeficientes seja computacionalmente eficaz, 
reduzindo o tempo gasto no processo. Vale ressaltar que a estratégia utilizada não 


foi vista em nenhum outro texto da literatura. 


Comparando os resultados obtidos com aproximações para o mancal curto 
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apresentadas porlFRISWELL et al./(2010) e utilizando a mesma simulação da Tabela 


obtém-se os seguintes erros relativos, tendo a aproximação como referência: 


Em Em o pd Pr] 


E. Fes [0.005 0.04 
Ex, Ex 10.001 0.04 


Ec. Ee |0.06 he EA 


Cyx Cyy 


Também é possível perceber, através das Figuras e 4.12] que os coefici- 





entes de rigidez e amortecimento obtidos pelo método apresentado são verificados, 
quando comparados as aproximações para um mancal curto, considerando diferentes 
cargas externas W aplicadas no eixo. Nos gráficos das Figuras [4.11]e [4.12] as linhas 
contínuas indicam os coeficientes obtidos através das aproximações apresentadas 


por |FRISWELL et al.| (2010), enquanto os termos obtidos pelo presente trabalho 


são indicados pelos símbolos apontados na legenda. 


1e7 Coeficientes de rigidez para diferentes valores de W 





0 100 200 400 500 


Wi NJ 


Figura 4.11: Coeficientes de rigidez para diferentes cargas W no mancal curto. 


Além da verificação do modelo proposto com a literatura, os gráficos das 
Figuras e apontam algumas outras informações. No que diz respeito aos 
coeficientes de rigidez da Figura4.11| é possível perceber que os termos diretos (kyy € 








kyy) aumentam para cargas externas maiores aplicadas, todavia k,, manifesta maior 
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1e6 Coeficientes de amortecimento para diferentes valores de W 





0 100 200 300 400 500 


WIN] 





Figura 4.12: Coeficientes de amortecimento para diferentes cargas W no mancal 
curto. 


variação a partir de W > 200 [N] do que kry. Já os termos cruzados (kr, € kyz) 


decrescem nessas mesmas condições. O coeficiente kz, apresenta pouca variação, 





entretanto k,y decresce rapidamente. A maior variação dos coeficientes de rigidez 
que envolvem a componente vertical das forças resulta do fato desta se equilibrar 


com o carregamento externo aplicado W. 


Para os coeficientes de amortecimento, o gráfico da Figura |4.12] mostra os 
termos diretos (Cry € Cyy) crescendo a medida que W aumenta. Mais uma vez, o 
coeficiente relativo à força vertical (c,,) é mais impactado do que o relativo à força 


horizontal (crr). Contudo, os termos cruzados decrescem de maneira semelhante, 


concordando com a afirmativa delFRISWELL et al.| (2010), que considera c;, = Cya- 


Os resultados preliminares obtidos neste capítulo demonstram que o modelo 
proposto se adequa, de maneira convincente, ao contexto de mancais hidrodinâmi- 
cos. Uma vez que é possível alcançar, através de processos simples, os principais 


componentes utilizados para análises da estabilidade de sistemas dinâmicos rotati- 
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vos, a modelagem revela-se satisfatória para o contexto supracitado. Contudo, vale 
destacar que a parametrização da geometria do problema, apresentada no Capí- 
tulo |2l deste texto, permite ampliar sua aplicabilidade a diferentes configurações de 


mancais, ou até mesmo para outros elementos de máquina. 


84 


5 RESULTADOS E DISCUSSÕES 


Este capítulo apresentará os principais resultados obtidos a partir de si- 
mulações computacionais para a modelagem proposta neste trabalho. O algo- 
rítmo utilizado foi desenvolvido em Python, sendo um dos resultados desta pes- 
quisa, e encontra-se disponível na plataforma GitHub, através do software Ross- 
Rotordynamicdl] uma biblioteca de código aberto para análise rotodinâmica, por 
(2020), no qual a autora deste texto possui coautoria. Atualmente 
neste software, o presente modelo e todas as suas funcionalidades correspondem à 
parte denominada Fluid Flow, que é responsável pela modelagem do filme de fluido 
em mancais e fornece ao resto do programa informações necessárias para o estudo da 
estabilidade de sistemas dinâmicos rotativos. Como Ross é uma biblioteca aberta, 
sendo constantemente atualizada, o link de acesso ao código utilizado neste texto 


encontra-se registrado também no Apêndice [4] 


Em razão do baixo custo computacional da solução numérica implementada, 
será viável avaliar os impactos da variação dos principais parâmetros que influenciam 
nos resultados do modelo, gerando uma ferramenta de análise para a construção 
de designs de operação. A fim de ilustrar algumas das possibilidades oferecidas 
pela descrição da geometria escolhida, serão gerados resultados para três diferentes 
configurações de mancais hidrodinâmicos: cilíndrico, elíptico e com desgaste. Nesse 


contexto, as alterações no campo de pressão, forças hidrodinâmicas, posição de 





equilíbrio e coeficientes de rigidez e amortecimento serão discutidas e interpretadas 


para cada caso. 


Todas os resultados que serão apresentados foram obtidos através de simula- 


https://github.com/ross-rotordynamics/ 


SO 


ções realizadas em um notebook com processador Intel Core 17 8565U Q 1.80GHaz, 
com memória RAM de 8GB Q 2400/2660 MHz e disco HD de 1 TB Q 5400 rpm. 


5.1 Mancal cilíndrico 


A estrutura cilíndrica de um mancal foi estudada no Capítulo |4] deste tra- 


balho. Sabe-se que a modelagem proposta se adapta de maneira satisfatória a este 





contexto, tornando possível realizar um estudo que objetiva entender o impacto da 


variação de alguns parâmetros nas respostas do modelo. 


Utilizando a simulação numérica da Tabela serão modificados, dentro 
de determinados intervalos, valores para a excentricidade, velocidade de rotação do 
eixo, carga externa aplicada, viscosidade do fluido e comprimento do mancal, para 
cada situação analisada. A partir dessas variações, será possível chegar a algumas 
conclusões acerca das condições geométricas e operacionais ideais para a estabilidade 
do sistema dinâmico rotativo, bem como sobre os fatores que mais impactam na 


resposta do modelo. 





Elemento Valor Variação 
Raio do estator (R,) 0.015 [m| - 

Folga (F) 9-107º [m] - 
Excentricidade (e) - 01F <e<o09F ml] 
Comprimento (L) 0.02 [m| 0.125 < L<012[m| 
Carga externa (W) 50 [N| I0<W < 150/N] 

Viscosidade (pt) 5.449. 1072 [Pa.s] | 0.01 < u <0.21[Pa.s] 
Rotação (tw) 52.36 [rad/'s| 50 <w< 550[rad/s| 


Tabela 5.1: Simulação numérica para o mancal cilíndrico. 
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5.1.1 Excentricidade 


A excentricidade, isto é, a distância entre os centros do rotor e do estator em 


um mancal cilíndrico, é a principal responsável pela variação na espessura do filme 





lubrificante e, consequentemente, pela formação do campo de pressão. Para que a 
descrição da geometria seja coerente com a realidade, a excentricidade nunca deve 
exceder o tamanho da folga radial, caso contrário, o rotor seria descrito com partes 
fora do estator. Por esta razão, utilzando a simulação numérica descrita na Tabela 
9.) as investigações acerca deste parâmetro serão realizadas considerando-o como 
porcentagens da folga. A fim de avaliar apenas o impacto da excentricidade, o valor 


do ângulo de atitude foi fixado em 5 = 7/7. 


O gráfico apresentado na Figura 5.I|deixa claro a relação de que quanto maior 
a excentricidade, maior será a variação da pressão. Para uma excentricidade de 0.1 


da folga, o pico máximo de pressão é de 9 - 10º [Pa], enquanto para 0.8 da folga foi 





de 6- 10º [Pa]. Este resultado já era esperado uma vez que, aumentando a distância 
entre os centros do rotor e estator, o fluido sofre maior compressão na área de 
menor espessura do espaço anular. Como consequência, também as forças de reação 
do fluido são alteradas com a variação dessa distância. Isso é perceptível na 
Figura [5.2] na qual os valores das forças radial N e tangencial T', descritas na Seção 
[4.3] deste texto, são menores para pequenas excentricidades quando comparadas aos 


resultados para excentricidades maiores. 


No que diz respeito aos coeficientes de rigidez e amortecimento, apresentados 
na Seção deste trabalho, as Figuras (a-b) indicam que estes também são 


resultados impactados pelas variações da excentricidade. E possível perceber que, ao 








passo que o rotor torna-se concêntrico ao estator, os coeficientes tendem a aproximar- 
se uns dos outros e também apresentam uma considerável redução em seus valores, 


o que pode gerar problemas para a estabilidade do sistema dinâmico rotativo. De 


ST 


Pressão ao longo de theta 
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Figura 5.1: Comportamento da pressão no meio do mancal cilíndrico para diferentes 
excentricidades. 
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Figura 5.2: Forças do filme fluido para diferentes excentricidades. 
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Figura 5.3: Coeficientes dinâmicos para diferentes excentricidades. 


acordo com MACHADO) (2006), um mancal será mais estável quanto maior for a 


diferença, em uma dada direção, entre os coeficientes de rigidez cruzados (kr € kyz) 
e diretos (kry € kyy), sendo o termo direto maior do que o cruzado. O autor também 
indica que a principal função do amortecimento é amenizar a resposta vibratória do 


eixo, e afirma que a máquina rotativa pode entrar em colapso caso o amortecimento 





seja nulo. Entende-se então que valores pequenos para o amortecimento incidem em 


uma maior vibração do eixo. 





Destes resultados, conclui-se que a excentricidade é um fator de alto impacto 
para a resposta dinâmica do mancal, e que valores muito pequenos desse parâmetro, 


no contexto de mancais cilíndricos, não são desejados em sistemas rotodinâmicos. 


5.1.2 Rotação 


Uma vez entendido que a pouca excentricidade torna o mancal cilíndrico su- 
jeito às instabilidades do sistema dinâmico rotativo, analisar fatores que influenciam 


na posição de equilíbrio do rotor passa a ser fundamental. Vale lembrar que a ex- 
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Posição do centro do rotor 
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Razão de excentricidade para diferentes rotações 





S = 1.59 
o = 550 
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Razão de excentricidade (e/F) 
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(a) Centro do rotor. (b) Variação da excentricidade. 


Figura 5.4: Variação da posição de equilíbrio para diferentes rotações. 


centricidade é uma característica gerada pelas forças de reação do filme fluido, e que 


estas, por sua vez, originam-se da rotação do sistema. 


As Figuras [5.4 (a-b) mostram, respectivamente, a posição do centro do rotor 
e a excentricidade para valores da rotação no intervalo 50 < w < 550[rad/s). Foi 
utilizada a simulação da Tabela na qual é dada uma carga W = 50|N| e o 


método iterativo para encontrar a posição de equilíbrio, apresentado na Subseção 





deste texto, é empregado. Analisando os dois gráficos é possível perceber que, 
como esperado, a medida que a velocidade de rotação fica mais intensa, o centro 
do rotor desloca-se em direção ao centro do estator e, por consequência, o valor da 
excentricidade diminui. Esse comportamento é similar ao apresentado na Figura 
4.81 que descreve o caminho percorrido pelo centro do rotor em função do Número 
de Sommerfeld (4.8). Rotações mais baixas geram um Número de Sommerfeld me- 
nor, enquanto para valores mais elevados, este adimensional torna-se maior. Desta 
maneira, o centro do rotor transita por uma curva próxima a um semicírculo rumo 


ao centro do estator. 
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Figura 5.5: Comportamento da pressão ao longo de O para diferentes rotações. 
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Figura 5.6: Comportamento das forças radial e tangencial para diferentes rotações. 
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Ea Coeficientes de rigidez Coeficientes de amortecimento 
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(a) Coeficientes de rigidez. (b) Coeficientes de amortecimento. 


Figura 5.7: Coeficientes dinâmicos para diferentes rotações. 





O comportamento da pressão para essas diferentes rotações é representado na 
Figura Aqui duas características se destacam: a altura e largura dos picos em 
cada caso. Para rotações menores, a curva do campo de pressão mostra-se estreita e 
com valores elevados, enquanto para rotações mais altas, a curva torna-se mais larga 
e com valores da pressão significativamente mais baixos (valores para comparação 
quantitativa). Por outro lado, o gráfico da Figura indica que os valores da 
força radial decrescem, para essas mesmas situações, ao mesmo tempo que a força 


tangencial cresce. 


Outros resultados a serem observados são os coeficientes de rigidez e amor- 
tecimento das Figuras (a-b). Com o aumento da rotação, a diferença entre os 
termos kz € kry diminui, chegando até mesmo à situação em que o termo cruzado 
é maior do que o direto, o que piora a estabilidade do sistema dinâmico rotativo. 
No entanto, os coeficientes k,, € kyz permanecem bem distintos em si e sempre com 
o termo direto maior que o cruzado. Já o amortecimento manifesta maior sensi- 
bilidade à rotação. É possível perceber que todos os seus coeficientes apresentam 


uma grande redução para rotações mais altas. Estes fatos, em acréscimo dos re- 
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(a) Centro do rotor. (b) Variação da excentricidade. 


Figura 5.8: Variação da posição de equilíbrio para diferentes cargas. 


sultados da pressão e força, podem ser justificados pela redução da excentricidade 
característica para altas velocidades de rotação. Ao girar mais rapidamente, o eixo 
deloca-se para o centro do estator, fazendo com que o fluido seja menos comprimido 
no espaço anular. Os resultados da rigidez e amortecimento indicam uma possível 
instabilidade no sistema dinâmico rotativo e, portanto, pode-se concluir que o man- 
cal cilíndrico talvez não seja a melhor escolha para o caso de máquinas que operam 


em altas velocidades. 


5.1.3 Carga externa 


Outro fator que influencia diretamente na posição do rotor é a carga externa 
W aplicada. Para avaliar a influência desse parâmetro na resposta do modelo, 


utilizou-se valores da carga no intervalo 10 < W < 150 [N] na simulação numérica 


da Tabela [,1] 
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Figura 5.9: Comportamento da pressão ao longo de O para diferentes cargas. 


Em oposição ao que ocorre quando a rotação aumenta, ao utilizar valores 
maiores para a carga externa W, o centro do rotor, ilustrado na Figura se 
afasta da origem. Isto é, a excentricidade aumenta a medida que o eixo torna-se 
mais pesado, conforme a Figura[5.8b] O Número de Sommerfeld (4.8) é inversamente 


proporcional a carga e, por essa razão, torna-se menor para cargas mais altas. 


Também o comportamento da pressão e das forças são influenciados por essa 
variação. À Figura |5.9/ indica que a largura das curvas da pressão não são tão alte- 
radas como ocorreu na variação da rotação. Contudo, o pico torna-se mais alto para 
eixos mais pesados, assim como as forças, expostas na Figura [5.10l Intuitivamente, 


estes resultados são coerentes com o esperado, uma vez que é necessário forças mais 





intensas para sustentar quaisquer objetos com pesos maiores. 








No que se refere aos coeficientes dinâmicos, as Figuras (a-b) apontam 
valores que diferem-se entre si a medida que W aumenta. Os termos diretos cres- 


cem para cargas maiores, enquanto os termos cruzados decrescem. Tanto para os 
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Figura 5.10: Comportamento das forças radial e tangencial para diferentes cargas. 





coeficientes de rigidez quanto para os de amortecimento, os coeficientes relativos à 
força vertical apresentam maior variação. Estes resultados sugerem que os mancais 
cilíndricos suportam elevadas cargas externas aplicadas, sem comprometer a esta- 
bilidade da máquina. Contudo, assim como para altas rotações, pode apresentar 


instabilidade do sistema dinâmico rotativo para eixos leves. 


5.1.4 Viscosidade 


Pela hipótese de um fluido newtoniano, considera-se a viscosidade como a 
constante de proporcionalidade na razão entre a tensão de cisalhamento e a taxa de 
deformação. Além disso, após as simplificações realizadas nas equações de Navier- 
Stokes, a equação diferencial parcial obtida pela modelagem apresenta este 
parâmetro como uma constante sendo inversamente multiplicada em dois de seus 


coeficientes: C4 (2.43) e C5 (2.44). Por esta razão, espera-se que este termo possua 
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Figura 5.11: Coeficientes dinâmicos para diferentes cargas. 





influência nos resultados da pressão e, consequentemente, nas forças, posição de 





equilíbrio e coeficientes dinâmicos. 





A simulação numérica empregada também é descrita na Tabela [5.1] com uma 
variação da viscosidade no intervalo 0.01 < yu < 0.21 /Pa.s], dentro do qual encontra- 
se os valores da viscosidade dos principais óleos lubrificantes utilizados em mancais 


nas indústrias. A variação da posição de equilíbrio para esses diferentes valores estão 





representadas nas Figuras nas quais percebe-se que óleos mais viscosos fazem 


com que o centro do rotor aproxime-se da origem, diminuindo a excentricidade e 


aumentando o Número de Sommerfeld (4.8). 





Consequentemente, também o campo de pressão e as forças do filme fluido são 
alterados por esse parâmetro, como pode ser visto nas Figuras Ç Fluidos 
menos viscosos geram um campo de pressão com valores mais elevados, entretando 
com uma curva mais estreita. Por outro lado, ao aumentar a viscosidade, a curva 
da pressão torna-se mais larga e com valores menores. Para as forças, a compo- 
nente tangencial cresce enquanto a radial decresce, ambas com um comportamento 


semelhante ao que ocorreu na variação da rotação (Figura |5.6). 
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(a) Centro do rotor. (b) Variação da excentricidade. 


Figura 5.12: Variação da posição de equilíbrio para diferentes viscosidades. 
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Figura 5.13: Comportamento da pressão ao longo de O para diferentes viscosidades. 
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Figura 5.14: Comportamento das forças radial e tangencial para diferentes viscosi- 
dades. 
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(a) Coeficientes de rigidez. (b) Coeficientes de amortecimento. 


Figura 5.15: Coeficientes dinâmicos para diferentes viscosidades. 
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Os coeficientes de rigidez da Figura mostram valores distintos entre si 





para a maioria das viscosidades, exceto para valores de 0.05 < yu < 0.075 |Pa.s], nos 
quais krrz, kry € ky tornam-se próximos. A partir deste intervalo, o termo cruzado 


(kry) da direção x torna-se cada vez maior que o termo direto (kr) da mesma 





direção. Como visto anteriormente, este é um indicativo de instabilidade no sistema 
dinâmico rotativo. O termo cruzado k,,, apesar de sofrer alterações em seus valores, 
permanece negativo e menor que o termo direto k,, em toda a variação. Já kry é 
negativo para fluidos com a viscosidade mais baixa, entretanto cresce rapidamente, 
tornando-se positivo, a medida que esse parâmetro aumenta. Para os coeficientes 
de amortecimento, a Figura [5.15b] indica que os termos diretos se intensificam para 


viscosidades maiores, contudo, os coeficientes cruzados apresentam pouca variação 





e são negativos em todo o intervalo. A partir destes resultados é possível entender 
que fluidos muito viscosos podem gerar instabilidades no sistema dinâmico rotativo 


devido à rigidez, entretanto tornam o sistema mais amortecido. 


5.1.5 Comprimento 


Por fm, o último parâmetro que será analisado é o comprimento L do man- 
cal. Como já foi visto no Capítulo |4] deste trabalho, este é um argumento muito 
utilizado para as habituais simplificações da literatura, uma vez que comprimentos 


infinitamente curtos ou infinitamente longos justificam o negligenciamento de uma 





das parcelas da equação de Reynolds (1.1). Torna-se interessante então entender o 
impacto deste fator nas respostas do modelo. Para tanto será utilizada a simulação 
da Tabela variando o comprimento a partir de sua relação com o diâmetro do 


mancal, no qual L/D será considerada no intervalo de 1/8 até 4. 


De acordo com os resultados expostos nas Figuras quanto maior for o 


comprimento do mancal, mais próximo os centros do rotor e do estator estarão, ou 
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Figura 5.16: Variação da posição de equilíbrio para diferentes comprimentos. 


seja, o Número de Sommerfel (4.8) aumenta e há uma diminuição na excentricidade. 


Também é possível perceber, através da Figura que os valores da pressão ao 





longo da posição angular diminuem consideravelmente a medida que L se torna 
maior, sendo este último fato justificado pela redução da variação da espessura do 


espaço anular. 


Já as forças do filme lubrificante apresentam dois comportamentos distintos 


no gráfico da Figura/5.18] Enquanto a força tangencial tem um rápido crescimento, a 





força radial diminui. Entretanto, percebe-se pouca variação desses valores a medida 


que o comprimento se aproxima de L = 4. 





E possível perceber na Figura [5.19alque os termos direto (kr) e cruzado (kr) 
da direção x são cada vez mais próximos até L/D = 1. A partir de então, a diferença 
entre estes coeficientes passa a aumentar, com o termo cruzado tornando-se maior 


que o direto, o que piora a estabilidade do sistema dinâmico rotativo. Para L/D > 1, 





também percebe-se que o termo k,, apresenta valores próximos a kxx, entretanto é 
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Figura 5.17: Comportamento da pressão ao longo de 6 para diferentes comprimentos. 
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Figura 5.18: Comportamento das forças radial e tangencial para diferentes compri- 
mentos. 
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Figura 5.19: Coeficientes dinâmicos para diferentes comprimentos. 


sempre maior que k,«. Em relação aos coeficientes de amortecimento, observa-se 
na Figura que os termos diretos crescem de maneira semelhante, enquanto 
os cruzados apresentam valores negativos, também próximos entre si. Entende- 
se então que o comprimento também é um fator que pode gerar instabilidades no 


sistema dinâmico rotativo, principalmente para valores da razão L/D > 1. 


5.2 Mancal elíptico 


De acordo comIMACHADO| (2011), os mancais de geometria cilíndrica, para 


casos de alta rotação ou com pequenas cargas externas aplicadas, ficam sujeitos à 
vibrações ou instabilidades do sistema dinâmico rotativo, o que também foi veri- 
ficado com os resultados da Seção anterior. Isso ocorre porque, nestas situações, 


a máquina tende a operar centrada, diminuindo a variação na espessura do filme 





lubrificante e, consequentemente, a capacidade de gerar pressão hidrodinâmica do 
mancal. O autor aponta que diversas configurações de mancais foram desenvolvidas 


a fim de resolver o problema com base nesta característica, produzindo pré-cargas 
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Figura 5.20: Desenho esquemático de um mancal elíptico. 


Fonte: Elaborada pelo autor. 


que fazem o eixo trabalhar fora dos novos centros estabelecidos, criando assim o 


efeito de cunha desejado. 


Segundo |OSTAY EN: BEEK! (2009), um componente comumente encontrado 


em máquinas de alta velocidade é o mancal elíptico ou 'mancal-limão”, como tam- 


bém é conhecido. Este é uma variação do mancal cilíndrico com ranhura axial e 


folga reduzida em uma direção. IOSTAY EN; BEEK! (2009) ainda afirmam que este 


modelo é relativamente barato e fácil de fabricar, pois é feito a partir de adapta- 





ções no próprio mancal cilíndrico, sendo formado por dois arcos circulares (cortes da 


estrutura original), com centros alinhados, gerando duas fendas horizontais, como 


indica a Figura 


Para a inclusão desta nova geometria, são necessárias adaptações no raio do 
estator, pois este não será mais constante em 6. Como observado na Figura O 
novo estator é composto pelo arco €4, com centro em O, unido ao arco (5, centrado 
em Os, ambos de raio R,. Nesta nova configuração, os centros estão a uma distância 


ce da origem, denominada elipticidade. E necessário descrever o estator a partir da 


103 





(a) m =0 (b)m5>1 


Figura 5.21: Variações do mancal elíptico de acordo com a pré-carga. 


Fonte: Elaborado pelo autor. 


origem. Esta nova distância será denominada Rj e apresenta variação ao longo da 


posição angular. 


Utilizando a lei dos cossenos, obtém-se a seguinte relação: 


Ri=vR?- esinfa+ecosa, (Da) 


(9) (9) 


T/2+0, se 0 € 1ºquadrante 
37/2+6, se 0 € 2ºquadrante 
0-T/2, se 0 € 3º quadrante. 
57/2-— 60, se 6 € 4ºquadrante 


onde a = 





Outro importante parâmetro a ser definido é a pré-carga m que, neste texto, 


será estabelecida como: 


(5.9) 


mM = 


Fº 





onde € é a elipticidade e E = R, — R, é a folga radial. 


Segundo |IMACHADO! (2006), a pré-carga pode ser considerada como a razão 


de não cilindricidade do mancal, variando entre valores maiores iguais a zero e me- 
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nores que um. Para m = O, o mancal torna-se cilíndrico, enquanto param > 1 os 


arcos do estator tendem a tocar no eixo (Figura |5.21). ID.R.GARNER; C.S.LEE; 
F.A.MARTIN! (1980) sugere que a pré-carga máxima adotada deve ser de 0.6. De 


acordo com os autores, valores acima de 0.7 podem ser utilizados no esforço de me- 
lhorar a estabilidade do sistema dinâmico rotativo, contudo, folgas laterais maiores 
são geradas, fazendo com que os coeficientes dinâmicos horizontais decresçam, além 


de exigir um maior fluxo de óleo. 


Elemento Valor 
0.015 [m| 

9: 107º [m] 
0.02 |m| 





5.449 - 1072 [Pa.s] 


| 
Carga (W) | 100 /N| 
) | 261.8 Jrad/s| 


Tabela 5.2: Simulação numérica para o mancal elíptico. 


Fonte: Adaptado de (MACHADO, |2006) 


Dada a simulação numérica da Tabela os resultados para a pressão no 
mancal elíptico, utilizando a modelagem apresentada neste trabalho, são compatíveis 
aos expostos na dissertação de (2006). As Figuras [5.22] (a-b) ilustram 
o campo de pressão para um mancal descrito como elíptico, porém com m = 0, ou 
seja, cilíndrico. Em contrapartida, os gráficos apresentados nas Figuras [5.23] (a-b) e 
5.24] (a-b), indicam o comportamento da pressão para pré-cargas iguais a 0.4 e 0.6, 


respectivamente. 


Observa-se que as Figuras C apresentam dois picos de pressão: um 
para o arco superior do mancal (0 a 7) e outro para o arco inferior (7 a 27). Este 
fenômeno pode ser explicado devido a configuração do mancal elíptico, pois, mesmo 


que o eixo opere centrado, haverá variações na espessura do fluido, fazendo com que 





o filme de óleo seja comprimido, tanto no arco superior, quanto no inferior. 
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Figura 5.24: Mancal elíptico com (m = 0.6). 
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Analisando também o gráfico da Figura entende-se que, quanto maior 


for a pré-carga aplicada ao mancal elítptico, maior será a força de reação gerada e, 





consequentemente, o campo de pressão. Este resultado também reforça a afirmação 


de |D.R.GARNER; C.S.LEE; F.A.MARTIN! (1980) pois, para m = 0.8, as pressões 


geradas são consideravelmente mais elevadas do que para as outras pré-cargas, sendo 





uma opção para casos em que se busca uma maior estabilidade dinâmica do rotor. 


te, Pressão ao longo de theta para o mancal elíptico 


Pressão [Pa] 





3 4 
Theta [rad] 


Figura 5.25: Pressão ao longo de 6 no meio do mancal elíptico para diferentes valores 
da pré-carga m. 


No que diz respeito à posição de equilíbrio, a Figura [5.26alilustra um compor- 
tamento diferente do que era visto no mancal cilíndrico, não sendo válida a relação 
do Número de Sommerfeld utilizada no mancal cilíndrico para analisar as caracterís- 
ticas de desempenho estático. Enquanto no primeiro caso, o centro do rotor traçava 
um caminho semelhante à um semicírculo, para o mancal elíptico variando a pré- 
carga, primeiro o rotor desloca-se verticalmente e, após m = 0.3, caminha na direção 


horizontal, aproximando-se do centro do estator, sempre reduzindo a excentricidade 


(Figura |5.26b). 
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(a) Centro do rotor. (b) Variação da excentricidade. 


Figura 5.26: Variação da posição de equilíbrio no mancal elíptico para diferentes 
pré-cargas. 
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Figura 5.27: Comportamento das forças radial e tangencial no mancal elíptico para 
diferentes pré-cargas. 
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Figura 5.28: Coeficientes dinâmicos do mancal elíptico para diferentes pré-cargas. 


As forças do filme fluido podem ser visualizadas na Figura |5.27) na qual é 


evidente que a força tangencial apresenta valores sempre muito maiores que a força 





radial para qualquer pré-carga aplicada no mancal elíptico. Também é perceptível, 
ao visualizar os gráficos das Figuras (a-b), que os coeficientes dinâmicos do 
mancal elíptico sofrem maior variação entre seus termos para valores maiores de m. 
Para a pré-carga m até 0.2, os termos kz, kry € kyy SÃO muito próximos, tornando a 
diferença entre eles imperceptível de visualizar no gráfico da Figura 5.2Yh. A partir 
de m = 0.3 é notável que a diferença entre os termos de uma mesma direção passa 


a aumentar, contudo o termo cruzado k,.,, é cada vez maior que o termo direto kg. 


xy: 


sendo este um indicativo de instabilidade do sistema dinâmico rotativo apontado 


por MACHADO! (2006). Já em relação aos coeficientes de amortecimento, todos os 


termos distanciam-se de zero a medida que a pré-carga aumenta. Vale lembrar que 


a principal função do amortecimento é amenizar os efeitos da vibração no mancal e 





que valores pequenos desses coeficientes não são desejados. 


Aqui fica claro a diferença entre os resultados do mancal cilíndrico e elíptico. 


No primeiro, todos os fatores que reduzem a excentricidade, implicam na diminuição 
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do campo de pressão, forças e variação dos coeficientes dinâmicos. Para o elíptico, 
ocorre exatamente o contrário. Ao variar a pré-carga, rotor e estator tornam-se 


mais próximos da concentricidade e, mesmo assim, isto significa aumento em to- 





dos os resultados citados. Entretanto, os coeficientes de rigidez mostram que este 
fato não necessariamente tornará o mancal mais estável, havendo ainda o risco de 


instabilidades no sistema dinâmico rotativo. 


5.3  Mancal com desgaste 


Apesar da lubrificação reduzir o atrito entre as superfícies metálicas do man- 


cal, estas estruturas costumam sofrer desgastes após um longo período operacional 





ou então devido a uma determinada quantidade de repetições dos ciclos de partida, 
onde existe maior atrito entre as partes sólidas por ainda não haver forças do fluido 
atuando no rotor. afirmam que uma das causas 
mais comuns de interrupções em sistemas rotativos é a ocorrência de falhas relaci- 
onadas a mancais hidrodinâmicos. Como o desgaste sofrido por estas componentes 


afetam a folga radial, certamente influenciam nos resultados do campo de pressão 





e, consequentemente, nas forças de reação e coeficientes dinâmicos. Por essa ra- 
zão, é relevante investigar o efeito do desgaste, em certas condições de operação, na 


resposta dinâmica do sistema. 


A geometria para o desgaste que será utilizada neste texto foi adaptada da 


versão apresentada poriIMACHADO; CAVALCA!|(2015). A descrição foi inicialmente 
proposta por DUFRANE; KANNEL; MCCLOSKEY| (1983) e validada experimen- 
talmente por HASHIMOTO; WADA; NOJIMA! (2008). 


DUFRANE; KANNEL; MCCLOSKEY!| (1983) consideram a espessura do 


desgaste uniforme na direção axial e localizada na parte de baixo do mancal, si- 
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Figura 5.29: Desenho esquemático de um mancal com desgaste. 


Fonte: Elaborada pelo autor. 


metricamente em torno do eixo Y. Porém, MACHADO; CAVALCA! (2015) afirma 


que, para casos práticos, a maioria dos desgastes ocorre no quarto quadrante, local 





da posição de equilíbrio do eixo para uma rotação no sentido anti-horário. Por esta 
razão, levando em conta que a região gasta é simétrica em relação a maior profun- 
didade do, será considerado um deslocamento angular 7 a partir do eixo vertical até 


a parte mais profunda, como indicado na Figura 


Para incluir o desgaste na geometria, é necessário realizar algumas adaptações 
no raio do estator. Considerando que a falha inicia na posição angular 6 = 0, e 


termina em 6 = 0, define-se a descrição do estator a partir da origem como: 


Ri=R,+do, (5.3) 


Õ, se 0<0<0, 0;<0<27m 
onde do — 
do—F(1I-+cos(0—-7/2)), se 0,<0<6; 


111 


Em 0, e em 0, a profundidade do desgaste é nula, então a localização das 


bordas pode ser definida da seguinte maneira: 


0, =— m/2+4 cos”! (do/F —- D) +79, 





0, =71/2-— cos! (do/F — 1) +49. (5.4) 
Elemento Valor 
Raio do estator (R5) | 0.015 [ml] 

Folga (F) |9-107º [ml] 
Comprimento (L) | 0.02 |m| 
Carga (W) | 18.9/N| 

Viscosidade (ju) | 0.1044 |Pa.s| 
Rotação (w) | 104.72 [rad/s| 


Tabela 5.3: Simulação numérica para o mancal com desgaste. 


Fonte: Adaptado de (MACHADO; CAVALCA! 2015) 


Com a geometria descrita, utilizando a modelagem deste trabalho aplicada 


aos dados da Tabela |5.3] os resultados do campo de pressão no mancal com desgaste 


são verificados com base nos estudos deIMACHADO; CAVALCA! (2015). Para um 


mancal cilíndrico sem desgaste, o campo de pressão é representado pelas Figuras 





(a-b). As Figurasf.31]|(a-b) ilustram o campo de pressão para o mancal com desgaste 
de profundidade máxima do = 10[um]) e um deslocamento y = 10º = 0.17[rad|. Nas 
Figurasf5.32](a-b) é considerado do = 50[um] para o mesmo deslocamento e, por fim, 
as Figuras [5.33] (a-b) mostram a pressão para o mesmo desgaste com o deslocamento 


y = 30º = 0.52frad). 





É possível perceber que, mesmo para pequenos desgastes, o campo de pressão 
é afetado drasticamente. Além de passar a apresentar dois picos, assim como ocorria 


no mancal elíptico, a magnitude da pressão também é alterada. Quando maior é a 
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Pressão ao longo de theta no meio do mancal. 






+ 
6 & 
60000 
o 
b., 
4 O 50000 
“o 
u 
3 Ú e 40000 
à o 
2 15 
o 30000 
1 à 
0 20000 
0.00 
“a cd 5 28 10000 
Om. 1.00 ERIC 
pj. 125 > 2.999 
Nena 75 1 ação ! 
2nj 2:00 0 pos 0 1 2 3 4 5 6 
Theta [rad] 


(a) Distribuição de pressão. (b) Vista lateral 


Figura 5.30: Mancal sem desgaste. 


Pressão ao longo de theta no meio do mancal. 
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(a) Distribuição de pressão. 


Figura 5.31: Mancal com desgaste (do = 10[um] e y = 10º). 


113 


Pressão ao longo de theta no meio do mancal. 
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(a) Distribuição de pressão. (b) Vista lateral. 


Figura 5.32: Mancal com desgaste (do = 50[um] e y = 10º). 


Pressão ao longo de theta no meio do mancal. 
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(a) Distribuição de pressão. 


Figura 5.33: Mancal com desgaste (do = 50[um] e y = 30º). 
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profundidade dy, maiores são os valores para a pressão e também mais elevado se 


torna o a primeira curva, manifestada pelo surgimento do desgaste (Figura |5.34). 


Pressão ao longo de theta 
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Figura 5.34: Pressão ao longo de 6 no meio do mancal com desgaste para diferentes 
diferentes profundidades do (y = 10º). 


A profundidade do desgaste na estrutura também influencia na posição de 
equilíbrio do rotor, como é visto nas Figuras Quanto maior é a profundidade 
do, mais o centro do rotor afasta-se do centro do estator, aumentando a excentrici- 
dade. Além da direção oposta ao centro do estator, o caminho percorrido também 
difere no comportamento quando comparado ao mancal cilíndrico sem desgaste, não 
sendo mais válida a relação com o Número de Sommerfeld (4.8). Para maiores pro- 
fundidades do desgaste, o centro do rotor localiza-se mais abaixo, descrevendo um 


percurso próximo de uma reta. 


Mesmo apresentando um comportamento do campo de pressão parecido com 
o do mancal elíptico, as forças radiais e tangenciais para o mancal com desgaste 
desenvolvem um comportamento distinto. O gráfico da Figura |5.36] indica que a 


componente radial da força sofre um rápido aumento para profundidades maiores 
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(a) Centro do rotor. (b) Variação da excentricidade. 


Figura 5.35: Variação da posição de equilíbrio no mancal com desgaste para dife- 
rentes profundidades do (y = 10º). 


do desgaste, ao mesmo tempo que os valores da componente tangencial descrescem. 








Finalmente, no que se refere aos coeficientes dinâmicos, as Figuras (a-b) 
descrevem o comportamento da rigidez e do amortecimento para diferentes profun- 
didades do desgaste. Há uma intensa variação dos coeficientes de rigidez, com todos 
os termos distintos entre si e apenas os valores de k,, negativos. Para desgastes de 
profundidade máxima até do = 30[um|, o termo cruzado k,,, é maior que o termo 
direto k,», indicando instabilidade do sistema dinâmico rotativo. A partir desse 


valor, k»y cresce, distanciando-se do termo cruzado, mas volta a decrescer logo em 





seguida. Os coeficientes relativos às forças verticais não apresentam variações muito 
relevantes, sendo o termo direto sempre muito maior do que o cruzado. Já os coe- 
ficientes de amortecimento desenvolvem valores positivos em seus termos diretos, e 
negativos e semelhantes em seus coeficientes cruzados. Em ambos os casos, não são 


estabelecidos padrões de crescimento ou decrescimento para seus termos. 


Assim, pode-se que concluir que desgastes no mancal cilíndrico são variações 
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Forças do Filme Fluido 


Forças [N] 
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Figura 5.36: Comportamento das forças radial e tangencial no mancal com desgaste 
para diferentes profundidades do (y = 10º). 

















1e6 Coeficientes de rigidez Coeficientes de amortecimento 
ij º o 25000 
* 
1.0 * » » 
* + 
+ E ii + 20000 
e 
0.54 
) 15000 
e kxx » + A + ” á 
a RR kxy 10000 
x Ú 
Tv kyx 
-0.5 + kyy 5000 
-1.0 + 0 
-1.5 7 á v E : v , ic 
Y v 
10 20 30 40 50 60 70 80 10 20 30 40 50 60 70 80 
dO [um] do [um] 
(a) Coeficientes de rigidez. (b) Coeficientes de amortecimento. 








Figura 5.37: Coeficientes dinâmicos do mancal com desgaste para diferentes profun- 


didades do (y = 10º). 


11% 


da geometria que afetam diretamente na resposta dinâmica do sistema. Os coe- 
ficientes de rigidez acusam intervalos em que o risco de instabilidades no sistema 
rotativo são mais críticos, o que pode servir para indicar condições operacionais que 


levariam o rotor a entrar em colapso. 


5.4 Comparações entre as geometrias 


Grandes partes das máquinas rotativas atuais, operam em altas velocidades 


de rotação. Como visto na Subseção [5.1.2] tal velocidade influencia em toda a res- 





posta dinâmica do mancal, desde a posição de equilíbrio até os coeficientes de rigidez 
e amortecimento. Por essa razão, é fundamental compreender como as respostas do 
modelo se comportam para variações da rotação em cada uma das geometrias apre- 


sentadas. 


Figura 5.38: Desenho esquemático relacionando as três configurações de mancais 
apresentadas. 


Fonte: Elaborada pelo autor. 
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Elemento Valor 
0.015 [ml] 
9:10" [m] 
0.02 [m| 





Raio do estator (R, 
Folga (F 
Comprimento (L 


| 
Carga (W) | 100 |N| 


0.05449 [Pa.s] 


Rotação (w) | 50 <w< 550 /rad/s| 


Viscosidade (ju 
( 
Tabela 5.4: Simulação numérica para a comparação das geometrias do mancal. 


A fim de comparar os resultados, foi utilizada a mesma simulação, cujos 
dados estão na Tabela [5.4] adaptando somente as especificidades de cada geometria 
e variando a velocidade de rotação no intervalo 50 < w < 550[rad/s|. Para o 
mancal elíptico considerou-se a pré-carga m = 0.4 e, para o mancal com desgaste, 
a profundidade máxima do = 50[um| com um deslocamento angular y = 10º. O 
desenho esquemático que relaciona as três configurações de mancais está ilustrado na 


Figuralo.38, na qual as geometrias são indicadas pelas letras CU, E e D, representando 





mancais cilíndrico, elíptico e com desgaste, respectivamente. 
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(a) Centro do rotor. (b) Variação da excentricidade. 


Figura 5.39: Posição de equilíbrio para diferentes configurações de mancais, variando 
a rotação. 
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As Figuras 5.39] ilustram a variação da posição de equilíbrio para o mancal 
cilíndrico, elíptico e com destaque. Nos três casos, a excentricidade diminui com o 
aumento da rotação, entretando é possível observar que estas aproximações ocorrem 
em caminhos diferentes. Enquanto o mancal cilíndrico desloca-se em um percurso 
semelhante a um semicírculo, como na Figura |4.8] do deslocamento em função da va- 
riação do número de Sommerfeld, o mancal elíptico gera uma curva mais acentuada, 
na qual a maioria dos pontos está alinhada ao eixo horizontal. Em contrapartida, o 
mancal com desgaste apresenta excentricidades bem mais elevadas para todas as ro- 
tações. É interessante destacar a diferença entre as posições de equilíbrio do mancal 
com e sem desgaste. Apesar de terem a geometria bem próximas, o desgaste altera 
muito a posição do centro do rotor e a forma como ele se comporta para diferentes 


rotações. 


Pressão ao longo de theta 
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300000 
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100000 





6 8 
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Figura 5.40: Pressões ao longo de O no meio do mancal para diferentes configurações 
da geometria, variando a rotação. 


Como já foi visto, o comportamento da pressão é alterado para cada geome- 


tria. A Figura /5.40| demonstra essas alterações de maneira de facilitar a visualização. 





Percebe-se que, além nas novas curvas, o valor máximo da pressão também é modi- 
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ficado. O mancal cilíndrico apresenta o menor valor da pressão, enquanto o elíptico 
atinge a maior magnitude, mesmo sendo a configuração com menor excentricidade. 
Também é possível perceber que o desgaste, além de criar uma nova curva, faz 
com que os resultados da pressão sejam mais altos do que no mancal cilíndrico, 


modificando também o local onde a pressão atinge seu maior valor. 


Forças 


Forças [N] 


cilíndrico 


elíptico re 
desgaste 





100 200 300 400 500 
w [rad/s] 


Figura 5.41: Forças para diferentes configurações de mancais, variando a rotação. 
A linha contínua representa a componente radial e a linha tracejada, a componente 
tangencial. 


As componentes radiais e tangenciais das forças do filme fluido (Figura [5.41), 
bem como os campos de pressão, são distintos para os três casos apresentados. 
Nos mancais cilíndricos e elípticos, a força radial N apresenta um decrescimento 
acentuado ao elevar a rotação, chegando até mesmo a alcançar valores negativos 
para o caso do mancal elíptico. Para o desgaste, o decrescimento dessa componente 
é moderado, entretanto com valores muito mais altos do que a força tangencial 
T, em crescimento também suave. A força tangencial 7 para o mancal elíptico 


apresenta pouca influência com a rotação, ilustrando um rápido crescimento apenas 
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nos primeiros aumentos do intervalo analisado. Em contrapartida, demonstra sofrer 


alto impacto para a geometria cilíndrica. 


Por fim, observa-se nos gráficos das Figuras (a-d) e (a-d) as vari- 
ações de cada termo dos coeficientes de rigidez e amortecimento para as diferentes 
geometrias no intervalo da rotação adotado. Algumas características destacam-se 


nesses resultados e serão comentadas a seguir. 





Sobre os coeficientes de rigidez percebe-se que, no mancal elíptico, os termos 
kry € kyy apresentam descontinuidades, crescendo até w = 300[rad/s| e, logo em 
seguida, há uma drástica redução em seus valores, o que não ocorre nas outras 
duas geometrias expostas. O termo k,, manifesta pouca variação com a rotação 
para o mancal cilíndrico, entretando decresce para o mancal elíptico e desenvolve 
valores bem mais elevados para a geometria com desgaste. Os termos kyy, kyz € 
kyy apresentam comportamentos de crescimento e decrescimento semelhantes nos 
mancais cilíndrico e com desgaste, apresentando a maior diferança entre eles no 


termo cruzado k. 


A fim de analisar as condições de estabilidade dinâmica rotativa nos diferentes 
mancais, é necessário comparar os coeficientes direto e cruzado da rigidez para uma 
mesma direção. Tal comparação torna-se mais compreensível quando observados os 
gráficos das Figuras [5.43] (a-b), que ilustram as diferenças (kr — kzy) € (kyy — Kya) 
no intervalo de rotações aplicado. Se estas diferenças são negativas, significam que 
o termo cruzado possui maior valor que o termo direto, indicando instabilidade no 


sistema dinâmico rotativo. O desejado é que estes valores sejam altos e positivos. 





Dessa maneira, é possível perceber que as diferenças, para as três configurações de 
mancais na direção x tentem a diminuir, entretanto o mancal elíptico apresenta um 
comportamento distinto dos demais, com o termo cruzado maior que o direto para 


rotações no intervalo 100 < w < 300 /rad/s). Já a diferença dos termos na direção 
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Figura 5.42: Coeficientes de rigidez para diferentes configurações de mancais, vari- 
ando a rotação. 
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16 Diferença entre os termos direto e cruzado na direção x 1Diferença entre os termos direto e cruzado na direção y 
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Figura 5.43: Diferença entre os termos direto e cruzado da rigidez em cada direção. 


y é positiva para toda a variação nas três geometrias, mostrando-se bem menos 


uniforme também no mancal elíptico. 


Diferente do que acontece na rigidez, o amortecimento não se mostra tão 
sensível à presença de falhas. Todos os coeficientes desenvolvem valores próximos 
para mancais cilíndricos e com desgaste, exceto pelos termos diretos cy, distintos 


para rotações mais baixas, mas que acabam aproximando-se ao passo que a rotação 





aumenta. A geometria elíptica resulta em coeficientes de amortecimento que seguem 
os mesmos padrões de crescimento e decrescimento dos demais, entretanto sempre 
com valores mais elevados. Os termos cruzados c;, € Cyy apresentam valores crescen- 
tes e negativos para todos os casos. Contudo, com rotações mais altas, tornam-se 
positivos apenas no mancal elíptico. Já os termos diretos decrescem em todos os 


casos sem atingir resultados negativos. 


Os resultados obtidos deixam claro o impacto da variação do tipo de man- 
cal nas respostas do modelo. Os mesmos parâmetros de entrada foram utilizados, 
havendo alterações apenas nas descrição da geometria, e mesmo assim as respostas 


apresentaram grandes mudanças. O mancal elíptico, introduzido como uma alterna- 
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Figura 5.44: Coeficientes de amortecimento para diferentes configurações de man- 
cais, variando a rotação. 
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tiva para corrigir às instabilidades do sistema dinâmico rotativo geradas pelo mancal 
cilíndrico, mostrou-se pouco eficiente uma vez que os coeficientes de rigidez apon- 
tam características de instabilidade do sistema. Por outro lado, houve o pretendido 


aumento do campo de pressão mesmo em situações estáticas em que excentricidade 





diminui, juntamente com maiores valores dos coeficientes de amortecimento, cuja 





função é reduzir a vibração do eixo. Quanto ao mancal com desgaste, percebe-se 
que, por menor que seja sua alteração, a falha compromete toda a resposta do mo- 
delo. O campo de pressão é alterado drasticamente, assim como as forças e a rigidez. 
Por outro lado, a falha no mancal aparenta não influenciar tanto no amortecimento 


do sistema. 


Dessa maneira conclui-se que a descrição da geometria é fundamental para 
o sucesso da modelagem. Por menor que seja a variação, a forma do mancal in- 
fluencia nas respostas obtidas. Vale destacar que a parametrização escolhida para 
a descrição do problema permite que tais alterações sejam feitas de forma natural, 
sem a necessidade de modificar o modelo apresentado, sendo possível ainda incluir 
novas mudanças que podem ser feitas através de usinagem ou observadas no próprio 
processo. Os resultados deste capítulo não classificam as configurações de mancais 
como eficientes ou não, mas sim constituem uma ferramenta para a escolha de me- 
lhores parâmetros, a fim de otimizar o funcionamento das máquinas rotativas dentro 
de seus contextos específicos. Utilizando o código FluidFlow as escolhas desses pa- 
râmetros podem ser feitas através da análise do desempenho desejado para o campo 


de pressão, forças, posição de equilíbrio e coeficientes dinâmicos. 
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6 CONCLUSÕES 





Considerando o escoamento de uma pequena espessura de filme de óleo lu- 
brificante entre dois cilindros, o presente trabalho apresenta uma revisão de toda a 
Teoria de Lubrificação, partindo das equações de Navier-Stokes até chegar a 
uma equação diferencial parcial elíptica que descreve o campo de pressão no 
espaço anular, correspondente à conhecida equação de Reynolds normalmente 


utilizada na literatura sobre o tema. 


Através de uma análise adimensional, foi possível promover um mapeamento 
dos termos das equações, ordenando-os por relevância no fenômeno físico para esse 
tipo de situação. Este mapeamento não é explorado neste texto pois segue-se as 


simplificações realizadas pela Teoria de Lubrificação clássica, entretanto fornece uma 





abertura para trabalhos futuros ao explicitar termos que possivelmente impactam 


mais na equação do que outros. 


Ao revisitar esta formulação, também foi possível escolher uma nova e eficaz 
parametrização da geometria, que viabiliza uma descrição do problema mais deta- 
lhada, permitindo adaptar a modelagem deste trabalho a outros contextos. Sem 
alterar o modelo proposto, realizou-se simulações em diferentes configurações de 
mancais hidrodinâmicos. Esta parametrização ainda permite que modelo seja adap- 
tado futuramente para muitas outras geometrias, incluindo alterações não somente 
na direção tangencial, como também na axial. Além disso, seguindo a literatura base 
C (2008). as equações foram trabalhadas 
em coordenadas cilíndricas com a finalidade de preservar os efeitos da curvatura, 
o que pode ser interessante para trabalhos futuros. Desta maneira, torna-se viável 


até mesmo realizar adaptações do modelo para outras componentes de máquinas 
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rotativas. 


Após as simplificações realizadas, o campo de pressão foi obtido numerica- 


mente através do método de Diferenças Finitas centradas, seguido de um adequado 


tratamento à esparcidade do problema, que resultou na redução do custo computa- 





cional envolvido no processo de resolução das equações que descrevem este campo. 
Tal resultado foi verificado a partir de sua aplicação em simulações para mancais 
hidrodinâmicos cilíndricos, comparando as respostas deste modelo com as obtidas 
através de aproximações clássicas da literatura, utilizadas na maioria dos trabalhos 


sobre o tema. 


Dentro da aplicação do modelo ao contexto de mancais hidrodinâmicos, ainda 
foi viável a obtenção de outras respostas do sistema, como as forças de reação do 
filme fluido, posição de equilíbrio do eixo e coeficientes dinâmicos. Como o campo 
de pressão é discreto, encontrou-se as forças a partir de um processo de integra- 
ção numérica, utilizando a regra de Simpson composta. Também foi proposto um 
método iterativo para a posição de equilíbrio que, através de Mínimos Quadrados, 
minimiza uma função residual para comparar a carga resultante com a carga externa 
aplicada. Por fim, foi apresentada um nova metodologia, diferente da comumente 
utilizada na literatura, para gerar os coeficientes de rigidez e amortecimento. Essa 
alternativa, através de uma solução também por Mínimos Quadrados obtida pelas 


pseudo inversas de Moore-Penrose, permite avaliar o impacto da perturbação ape- 





nas no momento em que se deseja calcular os coeficientes dinâmicos, tornando o 


processo computacionalmente eficaz. 


Outro relevante resultado é a construção do algorítmo em Python, desen- 
volvido pela autora desta pesquisa e utilizado em todas as simulações apresenta- 
das. Este código, denominado Fluid Flow, encontra-se disponível na plataforma 


GitHub, compondo a biblioteca para análise rotodinâmica de código aberto, Ross- 
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Rotordynamicg!] por TIMBÓ et al. (2020), no qual a autora deste texto possui 


coautoria. Através do Fluid Flow é possível escolher parâmetros geométricos e ope- 
racionais de acordo com análises do desempenho desejado da pressão, forças, posição 
de equilíbrio, coeficientes dinâmicos etc. Além disso, o baixo custo computacional 
da modelagem implementada permite gerar diferentes simulações, o que viabiliza 
sua utilização em algoritmos de aprendizado de máquina, função já explorada em 
Ross, que contém um módulo estocático para fornecimento desse tipo de análise. 
Este módulo também já explora incertezas acerca do modelo apresentado neste tra- 
balho como, por exemplo, o impacto da perturbação na posição de equilíbrio para 


a obtenção dos coeficientes dinâmicos. 


As facilidades geradas pelo baixo custo computacional deste modelo possi- 
bilitaram produzir análises sobre a variação dos principais parâmetros geométricos 
e operacionais do problema, expostas como resultados da pesquisa. Foram simula- 
das, em mancais cilíndricos, variações para a excentricidade, velocidade de rotação 
do eixo, carga externa aplicada, viscosidade do fluido e comprimento do mancal. 
Além disso, também apresentou-se outras duas configurações de mancais hidrodi- 
nâmicos: o mancal elíptico e com desgaste. Essas novas estruturas, obtidas através 
de mudanças apenas na descrição geométrica do problema, sem alterar a modela- 
gem, retratam o ganho resultante da parametrização escolhida. Em cada situação 
analisou-se o impacto dessas mudanças na resposta dinâmica do sistema e, por fim, 
foram feitas comparações entre as três geometrias, explorando as alterações nas 
respostas para diferentes velocidades de rotação. A diversidade das análises apre- 
sentadas constituem uma ferramenta que pode auxiliar em uma melhor escolha dos 
parâmetros, dentre as opções exploradas, para otimizar o desempenho da máquina 
rotativa. Contudo as análises realizadas têm caráter mais expositivo do que in- 
terpretativo, sendo possível ainda explorar mais formalmente o significado dessas 


variações nas respostas do mancal em relação à estabilidade do rotor. 


https://github.com/ross-rotordynamics/ 
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Em função de todos os resultados alcançados pela modelagem proposta neste 
texto, pode-se concluir que houve sucesso na abordagem. Vale ressaltar que as res- 


postas obtidas foram comparadas e verificadas com os principais métodos utilizados 





pela literatura sobre o tema. A presente pesquisa é relevante principalmente por 
sua característica de completude, uma vez que não é trivial encontrar trabalhos que 


sintetizem todos os assuntos tratados neste texto. Entretanto, ainda existem várias 





questões em aberto e considerações a serem feitas, indicando a oportunidade de 


trabalhos futuros. 
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APÊNDICE A CÓDIGO FLUID FLOW 


O algorítmo desenvolvido nesta pesquisa e utilizado nas simulações apresen- 


tadas foi desenvolvido em Python e encontra-se disponível na plataforma GitHub, 





compondo o software Ross-Rotordynamic4l] Atualmente neste software, o presente 
modelo e todas as suas funcionalidades correspondem à parte denominada Fluid 
Flow. Como Ross é uma biblioteca aberta e constantemente atualizada, o código 


exato utilizado nas simulações apresentadas neste texto, para fins de registro fixo, 


pode ser encontrado neste ou acessado pelo QR Code da Figura [4.1] 





Figura A.l: QR Code com link de acesso ao Fluid Flow: código exato utilizado para 
as simulações deste texto. 










https://github.com/ross-rotordynamics/ 


“nttps://github.com/ross-rotordynamics/ross/tree/8ce9846468fb86101eb2793dd30560 


7c29ca1l592/ 


